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CAPABILITIES AND APPLICATIONS OF THE 
PROGRAM TO OPTIMIZE SIMULATED TRAJECTORIES 
(POST) PROGRAM SUMMARY DOCUMENT 

By G. L. Brauer, D. E. Comick, 
and R. Stevenson 

Martin Marietta Corporation 
P. 0. Box 179 
Denver, Colorado 80201 


SUMMARY 

This report summarizes the capabilities and applications of the three- 
degree-of-freedom (3D0F) version and the six-degree-of-freedom (6D0F) version 
of the Program to Optimize Simulated Trajectories (POST). The document sup- 
plements the detailed program manuals (ref. 1, 2, and 3) by providing additional 
information that motivates and clarifies basic capabilities, input procedures, 
applications and computer requirements of these programs. The information will 
enable prospective users to evaluate the programs, and to determine if they are 
applicable to their problems. 

The report is presented in five chapters with the first containing a sum- 
mary of what is believed to be the important features of both programs. Chapter 
1 is intended to provide enough information to enable managerial personnel to 
understand the capabilities of the programs. The remaining chapters are pro- 
vided to describe the POST structure, formulation, input and output procedures, 
sample cases, and computer requirements. These chapters provide answers to 
basic questions concerning planet and vehicle modeling, simulation accuracy, 
optimization capabilities, and general input rules. Several sample cases are 
also presented. These sample cases contain enough detail to enable them to 
serve as guidelines for new users. Should more detailed questions arise, it 
is recommended that the POST Formulation Manual, Utilization Manual, and Pro- 
grammer's Manual be consulted. 


INTRODUCTION 

The original 3D version of POST was developed in 1970 as a Space Shuttle 
Trajectory Optimization Program. Since that time, the program has been signifi- 
cantly improved with additional capabilities added in the areas of vehicle 
modeling, trajectory simulation, and targeting and optimization. The program 
is capable of simulating and optimizing trajectories for a wide variety of 
aerospace vehicles operating in the vicinity of a single planetary body. 

The popularity of 3D POST led to the development of the 6D0F version. 6D 
POST is identical in design and use as its 3D0F counterpart. The principal 
feature of 6D POST, in comparison to other 6D0F programs, is it's easy to use 
input procedure. This capability was obtained by the development of a special 



NAMELIST input processor that does not contain the namelist size limitation of 
the standard NAMELIST routine. 

During the development process, considerable effort was made to ensure 
that both versions were generalized in capability, yet easy to learn and use. 

As a result, 3D and 6D POST can be readily used by trajectory engineers without 
specialized training in areas such as optimization theory. Ease of use and the 
ability to be used on almost any kind of near-Earth trajectory problem has re- 
sulted in considerable interest in these programs throughout the industry. This 
interest has resulted in numerous questions being asked by potential users con- 
cerning the general capabilities of both programs. Thus, the purpose of this 
summary report is to answer these kinds of questions without requiring reference 
to the detailed program manuals. As a result, this report contains a broad 
spectrum of information related to program features, structure and design, 
utilization, sample cases, and computer requirements. These data will provide 
the potential user with basic capability information, and the experienced user 
with a summary for quick reference purposes. 

The instruction manuals and source tapes for both 3D and 6D POST are avail- 
able from: 

Computer Software Management & Information Center 
112 Barrow Hall 
University of Georgia 
Athens, Georgia 


PROGRAM FEATURES 

The 3D POST and 6D POST are general purpose- FORTRAN codes designed for 
flexible 3D and 6D simulation and optimization of trajectories for aerospace 
vehicles. A summary of the key program features is presented in Figure 1. 

In reviewing the program features, it is important to realize that 3D POST 
and 6D POST are separate programs. However, the executive structure and I/O 
characteristics are identical. The only significant difference between these 
two programs is that the rotational equations of motion are included in 6D POST 
as depicted in Figure 2. It is important to note, however, that the 6D version 
requires an additional 44 131 octal storage locations over its 3D counterpart. 


Simulation Capabilities 

POST is best described as an N-phase trajectory simulation program. This 
means that the POST input processor and executive structure enables the user to 
simulate the trajectory by a logical sequence of trajectory phases. In each 
phase, physical and nonphysical aspects of the simulation can be modeled to any 
accuracy deemed appropriate by the trajectory analyst. In this manner, each 
phase of the trajectory can be simulated accurately and efficiently by appro- 
priate user input and program option selection. Figure 3 illustrates the rela- 
tionship between trajectory phases, events, and POST input data structure. In 
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PROGRAM FEATURES 

• Ea*y co l«arn SAMEL I ST input for 
both 3D and 6D siBulatlon • Auto- 
matic input error diagnostics • 
Hollerith specification of the 
targeting and optimisation formu- 
lation • Input trajectory by phase 

e Program operational on CDC 6000 
series, IBM 370, US (VAC 1108 and 
1110 • User-defined block print- 
out e Input echo e Optional phase 
summary print e optional print 
for special purpose calculations 
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OPTIMIZATION FEATURES 

• Discrete Parameter Targeting and 
Optimisation 

a Ability to select the optimiza- 
tion criteria, the conetralnts, 
and the control variables from 
a dictionary of over 400 program 
variables 

• Equality and inequality con- 
straints • Selection of several 
popular optimization algorithms e 
Automatic problem acal ing • Rapid 
calculation of numerical deriva- 
tives • Detailed and summary 
optimization search printout 
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SIMULATION FEATURES 

e Planet model selection inc ludes the 
1960 Fleher Earth model and the 
Smithsonian Earth model • Atmosphere 
mod«-l s.„*«?ction Includes the 1962 
U.s. Standard, the 1963 Patrick APB, 
or a general table model Input by 
the user • Winds e General weight 
and mass properties Input options • 
Rocket or alrbreethlng propulsion • 
All standard aerodynamic coefficient 
Input options * Laminar and turbu- 
lent aeroheatlng models • Five pop- 
ular attitude reference options • 
General GN6C moduals 
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TRAJECTORY FEATURES 

e Simulation by phase with 

unlimited number of trajectory 
phases • Generalized event 
stopping capability ♦ Primary 
secondary, roving, or repeat- 
ing events can be specified 
enabling any trajectory sequence 
of events to be modeled directly* 
Inlc lallzst ion trajectory and 
attitude in all popular reference 
systems 

• Several standard numerical 

integration methods * Keplerian 
and Encke orbital propae at 
options • Instantaneous ve loci tw- 
addle ions • Static t.rl\a pit:! 
and yaw • Vertical, rsko o::' hold- 
down model a Horizontal takeoff 
model 


Figure 1.- Summary of Program Features 


POST, the data structure is arranged according to the sequence of events defining 
the trajectory phases. It is noticed that every phase is specified by an event 
that defines the initiation of the phase. Thus, each event (other than the 
first and last events) serves two purposes: (1) it ends the previous phase, 

and (2) it starts the subsequent phase. The input data cards for each phase 
are located between the two events that define the phase. This basic input 
rule is also illustrated in Figure 3, where any or all of the simulated data 
can be changed in any phase. The data used in the simulation of a specific 
phase are the sum of the data input in all previous phases plus any new addi- 
tional data and/or options. If the added input data correspond exactly to any 
previous input data, then the program will use the latter data in the trajec- 
tory simulation. There are no restrictions on the number of events in a given 
problem, and the event criteria can be selected as any variable computed in the 
program. The capability to define the problem by phases is an important fea- 
ture of these programs because it enables complex problems to be formulated in 
a step-by-step fashion. 

The simulation capabilities can be categorized in three types: (1) the 

planet model, (2) the vehicle model, (3) trajectory simulation options. Sum- 
maries of these capabilities are presented in Tables 1, 2, and 3, respectively. 
Generally each particular simulation model has several options available. How- 
ever, if a particular model is not available automatically then general models 
can be used augmented by user input data. For example, two commonly used 
atmosphere models are available that require no user input. These are the 1962 
U.S. Standard and the 1963 Patrick AFB atmosphere models. However, in raanv 
cases, such as a Mars entry study, the user can and must input his own me. 
via generalized input of pressure, density, and temperature (or speed of gout.-' .• 
The detailed procedures for this type of input are described in reference 1. 
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Figure 2.- Relationship Between 3D and 6D POST 
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Figure 3.- Relationship Between Trajectory Sequence of Events 
and POST Input Data Structure 



TABLE 1.- SUMMARY OF PLANET SIMULATION OPTIONS 


Planet Model 
Characteristic 

Options 

Description 

Oblate Spheroid 

1. Spherical 

These models are defined by the equatorial radius, 

and 

2. 1960 Fisher 

the polar radius, the rotation rate, the gravita- 

Gravitation Model 

Earth 

tional constant y, and the second, third, and 
fourth gravitational harmonics, J 2 , J 3 , and J 4 . 


3. Smithsonian 
Earth 

For a spherical planet, J 2 , J 3 , and J 4 are all zero. 

Atmosphere 

1. General 

1. Atmospheric pressure, temperature* density, and 
speed of sound are computed from user input 
tables as a function of altitude. Speed of 
sound and density tables can be omitted, in 
which case they are computed from input values 
of ratio of specific heats, molecular weight* 
and the Universal gas constant. 


2. 1962 U-S. 

Standard 

2. The 1962 U.S. Standard Atmosphere model is 
given as a function of geopotential altitude. 
In POST, the molecular weight is assumed con- 
stant . 


3. 1963 Patrick AFB 

3. In this model, pressure and temperature are 
calculated as a function of geometric altitude 
and a set of prestored polynomial coefficients. 

Winds 

1. Geographic 

1. The wind velocity is input directly in the 
geographic frame by defining the Easterly, 
Northernly, and radial components of the wind 
velocity . 


2. Meteorologic 

2. The wind velocity is computed from the total 
wind speed and meteorologic wind heading. 


TABLE 2.- SUMMARY OF VEHICLE SIMULATION OPTIONS 


Vehicle 

Characteristic 

Options 

Description 

Propulsion 

1. Rocket Engines 

1 Vacuum thrust, nozzle exit area, and flowrate or 

*sp can k* Input for up to 15 engines per phase. 


2. Jet Engines 

2. The ratio oi total thrust to the atmospheric pres - 
sure ratio and the ratio of specific fuel consumptior 
to the square root of the atmospheric temperature 
ratio can be input for up to 15 engines per phase. 


3. Ramjet Engines 

3. Thrust coefficient and specific fuel consumption 
can be input for up to 15 engines per phase. 

Aerodynamic 


1. Lift, drag,- and sideforce coefficients (3D), and 
pitch, yaw, and roll aerodynamic coefficients 
(6D) are input as tables. 



2. Axial, normal, and sideforce coefficients (3D), 
and pitch, yaw, and roll aerodynamic coefficients 
(6D) are input as tables. 

Aeroheating 

1. Chapman’s Equation 

1. Heat rate is computed from Chapman's equation 
wi til t lie nose radius as an input. 


2. General Table 
Lookup 

2. Heat rate is computed as a table lookup based on 
as many as three Independent variables. 


3. Modified Chapman's 
Equation 

3. Heat rate is computed as the product of Chapman's 
heatrate and a general table lookup. 


4. Turbulent Flow 

4. Similar to modified Chapman's equations with 
different exponents. 


5. Maximum Centerline 

5. Heat rate is computed based on correction for 
altitude and angle of attack. 

Steering 

1 . Open Loop 

1. 3D attitude Is calculated from tables, polynomials, 
or linear feedback; 6D attitude commands are com- 
puted from tables. 


2. Closed Loop 

2. All closed loop models must be coded specifically 
for each application. 
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TABLE 3.- SUMMARY OF TRAJECTORY SIMULATION OPTIONS 


Simulation Faaturaa 

Option* 

Description* 

Integration Method* 

1. Rung* Kutta-4th Order 

1. The standard 4th order Runga-Kutta single atep 
method for a set of first order ordinary differ- 
ential equations. 


2. Rung* Kutta-8th Order 

2. An eighth order single step method. 


3. Predictor-Corrector 

3. A variable-step variable-order method developed 
by F. T. Krogh of the Jet Propulsion Laboratories. 

Table Interpolation 

1. Piecewise Constant 

1. The function is assumed to be a step function 
based on the user's input. 


2. Plec«wlse Linear 

2. Linear interpolation is used between data points. 
The tables can be monovariant, bivariant, or tri- 
variant. 


3. Piecewise Cubic 

3. Cubic interpolation is used between data points. 
The tables can be monovariant, bivariant, or tri- 
variant. 

Event* 

1. Primary 

1. tvents that must occur, and must occur in ascend- 
ing order according to the event number. 


2. Secondary 

2. Events that must occur in ascending order between 
their bounding primaries. The occurrence of a 
primary nullifies the previous secondaries. 


3. Roving 

3. Events that can occur any time after the occur- 
rence of all smaller primaries. 


4. Repeating 

4. Events that can repeat at a specified increment 
or at specified values. 

Orbital Propagators 

1. Laplace's Method 

1. A semianalytical method for propagating the posi- 
tion and velocity of a nonthrusting vehicle in 
vacuum flight over a spherical planet. 


2. Encke's Method 

2. A rapid method for propagating orbits perturbed 
by planet oblaceness and atmosphere. 

Special Purpose 

lT -V Additions 

1. At any specified event an Instantaneous velocity 
change can be added to the -vehicle. 


2. Static Trim 

2. The engine gimbal angles or flap deflection* sre 
computed to balance the pitch and yaw moments 
caused by thrust of other engines and aerody- 
namic forces. 

Launch 

1. Vertical Holddown 

1 , This model is used to simulate vehicle holddown 
by maintaining the position and velocity relative 
to the planet constant until launch. 


2 , Horizontal Takeoff 

2. This option is used to simplify horizontal take- 
off by allowing the vehicle to accelerate only 
in the local horizontal plane 


One additional degree of complexity usually arises in 6D trajectory work, 
that is, the guidance, navigation, and autopilot equations must generally be 
coded in FORTRAN and added to 6D POST in the appropriate subroutines. Proce- 
dures for accomplishing this effort are described in reference 2. 


Targeting and Optimization Capabilities 

POST has a complete discrete parameter nonlinear programming capability. 
This means that POST can minimize (or maximize) a user-selected performance 
function, subject to nonlinear target conditions and/or inequality constraints 
The control variables can be any parameter that influences the performance cri 
teria and/or the mission constraints. The performance function, the target 
conditions and constraints, and the control variables can be selected from a 
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dictionary of over 400 program variables. Some typical examples of these 
variables are shown in Figure 4. POST also contains several popular opti- 
mization algorithms that can be selected by the user. A brief summary of these 
algorithms is given in Table 4. As indicated, the accelerated projected gradi- 
ent algorithm is normally used as the basic optimization technique. This algo- 
rithm is a combination of Rosen's projected gradient method for nonlinear pro- 
gramming and Davidson's variable metric method for unconstrained optimization. 

In the targeting mode, the minimum norm algorithm is used to satisfy the tra- 
jectory constraints. The cost and constraint gradients required by these algo- 
rithms are computed normally as first differences calculated from perturbed 
trajectories. In some particularly difficult cases, symmetric differences are 
used to more accurately approximate the derivatives. To reduce the costs of 
calculating numerical sensitivities, only that portion of the trajectory in- 
fluenced by any particular independent variable is reintegrated on the perturbed 
runs. This feature saves a significant amount of computer time when targeting 
and optimization is performed. 



Typical Optimisation Variables 

Payload Weight, Burnout Weight, Launch Weight, 
Total Propellant Used, Burn Tines, 

Burnout Velocity, Burnout Altitude, 

Down rang*, Crosa range, Total lUnge. 

Typical Target Conditions and Conetralnta 

Altitude, Radius, Velocity, Flight Path Angle, 
Apogee/Perigee Altitude /Radius, 

Orbital Eleaents, Dovnrange /Cross range , 
qa, Total a, dynamic pressure, acceleration, 
Temperature, Keatrate, Total Heat Load. 


with respect to u 



Typical Control (Independent) Variable 

Attitude Angles (a, 6, o or 4>, 6, ♦), or Coefficients 
of Altitude Polynomials, Throttle Settings* 

Event Criteria Values (Burn Time, etc), 

Thrust Levels, Weights, Initial Conditions. 


Figure 4.- Optimization Summary 


TABLE 4.- TARGETING AND OPTIMIZATION PROBLEM SOLUTION METHODS 


General Class of Problems 

Examples 

Available Algorithm in POST 

Recommended Algorithm 

Constrained Optimization 
with Inequalities 

Optimize: f (u) 

Subject to: c_ (u) 0 

Ascent to Orbit 

Max: U 1M n (Payload V.’eigliL) 

Subject to: h “ 10*1 K05 ft 

(Altitude) 

'BO ' 0 < ri ’ A) 

V BO ’ » 

(Velocity) 

q • 400 psf (Dynamic 
Pressure) 

Accelerated Projected Grad- 
ient Method 

Single Penalty Functions 
Using: 

- Steepest Descent 

- Conjugate Gradients 

- Davidon’s Method 

Accelerated Projected 
Gradient Method 

Unconstrained Optimization 
Optimize: f(u) 

Sounding Rocket 

Max: h , .. (Max Altitude) 

h“0 

Steepest Descent 
Conjugate Gradients 
Davldon's Method 

Davldon's Method 

Target ing : 

Determine u Such That 
£ (u) - 0 

Space Shuttle Entrj/ 

Determine the F.ntrv y and Azimuth 
Such That Landing Site is Reached, 
i . e . 

Latitude ♦ - Specified 

Value 

Longitude q - Specified 
Value 

Steepest Descent 
Conjugate Gradients 
Davidon's Method 
Newton Raphson 

Newton Raphson 


8 


REPRODUCE TLITY OF THE 
ORIGINAL PAGE £3 POOR 






Input/Output Capabilities 


All program inputs are made using an improved NAMELIST input processor. 

The key features of this extended NAMELIST capability are illustrated in Figure 
5. The permissable sequence of NAMELISTs are also given in Figure 6. As indi- 
cated, $SEARCH and at least one $GENDAT are required per problem. $TBLMLT and 
$TAB are optional within $GENDAT depending on whether tabular data are required 
in that phase. 


REPEATED HOLLERITH SPECIFICATION 



SPECIAL NAMELIST LIST OPTIONS 

BLANK IN COLUMN 1 — LIST ONLY THOSE CARDS CONTAINING ERRORS 

P IN COLUMN 1 — LIST COMPLETE INPUT FILE 

L IN COLUMN 1 — LIST COMPLETE INPUT FILE INCLUDING A 

CARD COUNT 

Figure 5.- Extended Namelist Capabilities 
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Targeting & 

Optimization 

Input 


General Simulation 
Input for Phase 1 


Table Data 
for Phase 1 


General Data 
for Remaining 
Simulation 
Phases 


PI Sf arch 

C INPUT IN THIS NAMELIST ALL DATA PERTAINING TO S 

c 

C l) OPTIMIZATION AND TARGETING FORMULATION 

C optimization variable 

C CONSTRAINTS AND TARGET CONDITIONS 

C CONTROL VARIABLE 

C ?} ALGORITHM SELECTION AND CONTROL FLAGS 

ft 

PftGENDAT EVENT * 1, 

C INPUT IN TMIS NAMELIST ALL DATA AND PROGRAM CONTROL FLAGS PF AT A I NINC 
C TO: 

C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 


1> PLANET MODEL - DELATE SPHEROID 

- GRAVITY COEFFICIENTS 

- ATMOSPHERE MODEL 

2 ) VEHICLE MODEL- MASS PROPERTIES AND REFERENCE 
GEOMETRY 

- ENGINE OPTION ANG ENGINE LOCATIONS 

- AERO OPTION 

- GUIDANCF AND CONTROL 
2 TRAJECTORY - INITIAL CONDITIONS 

- INTEGRATION METHOD 

- SPECIAL OPTION CONTROL FLAGS 


PtTBLMLT 

C INPUT IN THIS NAMELIST ALL TABLE SCALE FACTORS FOR PHASE 1.0 

ft 

C INPUT THE THRUST TABLE C FOR PHASE 1.0 ) 

ft 

C INPUT THE NOZZLE EXIT APEA TABLE 

ft 

C INPUT THE CD OR THE CA TABLE FOR DROPTANK PLUS ORBITCR 

I 

C INPUT THE CL OR THE CN TABLE F OF DROPTANK PLUS ORBITIR 
C 

ENDRHS * 1, 

ft 

PftGENCAT EVENT * Z.O, CRITP * 4HASMG* VALUE *0.25, 

C TURN OFF HOLODOWN MOOE L VIA DATA INPUT HERE. 

ENDRHS * 1, 

ft 

RftGENDAT EVENT * 3.0 , C R I TR * 4HVFLR t VALUE 

ENDPHS *1, 

S 

RftGENDAT EVENT *A.O, CR I TR 

C INPUT DATA REQUIRE TO US E Cl 
ENDRHS xl f ^ "VaLUF « 6B00.0, 

PftGFNrAT EVE^I 

IHANGED SCALE FACTORS 

CD OR CA TABLES FOR ORBITER ( WITHOUT DROPTANK | 

C INPUT CL OR CN TABLES FOR OPE I TER I WITHOUT DROPTANK ) 

C 

ENDPHS *1, 

ft 

PftGENDAT EVENT =7.0, C R 1 TR = 5HTDVPP, VALUE = 5.0, 

C INRUT DATA REQUESTING LINES INF STEERING 
EDNPHS *1, 

ft 

PftGENDAT EVENT *6.0* CP I TR = AHA S MG « VALUE * 3.0, 

C INPUT DATA REQUESTING THROTTLING TO 3G LIMIT 
ENDPHS *1, 

ft 

PftGENDAT EVENT =R.O, CRITF * 5HWPR0P , VALUE « 0.0, 

C INPUT DATA TO TURN CF F ENGINES 
ENDPHS *1, 

ft 

PftGENDAT EVENT =10.0, C R 1 TP * 5HTDVPP, VALUE * 60., 

ENOPHS * 1 , 

FNDPRB =1, 

END JOB *!* 

ft 



Figure 6,~ POST Input Structure 


10 


The program also has the capability to accept input in either English or 
metric units. The only restriction being that the input units cannot be mixed. 
The output variables can also be printed in either English or metric units, 
but not both. 

All English to metric conversion constants can be changed by input if de- 
sired. The stored values of these constants were obtained from reference 3. 

Care must be taken when using metric system input and output. The input 
units must be of the same type as the English units. For example, values for 
weights must be input in units of force (newtons) rather than mass (kg) . Vari- 
ables output in nautical miles in English units must be kilometers in metric. 

The program also has the capability to print a summary of the table data 
input. This feature is useful in that input errors in the table format are 
easily detected by reviewing this printout. 

In addition to the normal input procedures described above, the program 
has the capability to process more than one input problem per pass at the com- 
puter. This feature, referred to as Multiple Runs, enables the user to change 
the basic input deck (which represents a single problem) by adding a subsequent 
set of data cards representing the changes (addition and/or deletions) imme- 
diately after the original set of data cards. Once the first problem is com- 
pleted, the program will automatically modify the input file as defined by the 
additional cards and run the resulting data file. This capability is extremely 
useful in performing parametric studies where only a few input variables are 
changed from one run to the next. 

There are two basic categories of POST output: (1) trajectory data, and 

(2) program control data. Trajectory data can be output in the form of a com- 
puter printout and/or as a profile tape. The typical computer printout con- 
sists of (1) an input echo, (2) an input summary, (3) a trajectory printout 
that optionally includes a phase data summary at the beginning of each event, 
(4) special trajectory printouts, such as orbital parameter, tracking data, 
etc, and (5) targeting and optimization iteration traces and summaries. Tra- 
jectory data can also be written on a profile tape for storage or auxiliary 
output purposes. Typical output obtained from the profile tape includes tra- 
jectory plots, orbital ground tracks superimposed on a world map, and punched 
cards. The computer routines required to generate these types of output from 
the profile tape are computer /system-dependent and not contained in POST, but 
are readily available at most modern computing centers. The second type of 
output, namely, program control data, is contained in the initial input summary 
and is updated in the phase data summaries output at the beginning of each 
event. Employing these summaries, the user can always determine exactly what 
options are being employed in the simulation during any phase of the trajectory 
simulation. The various forms of POST output are summarized in Figure 7. 
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POST Output Options 


Auxilary Options 
Using Profile Tape 


I 


II 


III 


Printed Output 
Options 


Profile 

Tape 

Option 


Graphics 

Output 



Input Echo 


Input Sumnarles 


Phase Data Summaries 


Block Printout 


Iteration Traces 


Iteration Sumnarles 



Plots 



or Cards 



Figure 7.- POST Output Options 


The basic trajectory printblock can be defined, in terms of size and con- 
tent, by the user. Any common variable computed in the program can be output 
by including its hollerith name in the printblock definition. If the user 
does not wish to define a printblock, then a nominal default printblock can be 
used. Any printblock may be modified with additions and/or deletions according 
to the rules described in reference 1. 


12 



The targeting and optimization output has two levels of detail. The 
iteration summary is the output most frequently used and it is always dis- 
played in the printout. This summary contains all the information required to 
understand the progress of the iteration process. Information, such as tra- 
jectory sensitivities, stepsize limits, univariant curvefit summaries, current 
control variables, current constraint errors, and current performance are in- 
cluded in the summary. When the iteration is judged to be improperly converg- 
ing, a detailed iteration trace can be requested. This trace gives all in- 
formation required to determine the exact iteration cycling calculations and 
prints the data in the exact chronological order that they were calculated 
internally. The data are usually only required when new problems are being 
formulated, and are extremely useful In identifying problem formulation and 
setup errors. It also gives key information that can be used by experienced 
POST users to speed up the convergence of the targeting and optimization algo- 
rithms . 


Program Applications 

During the past several years, POST has been used to solve hundreds of 
performance and mission analysis problems for atmospheric and orbital vehicles. 
A brief summary of the more standard applications of POST is contained in 
Table 5. The spectrum of problems shown in Table 5 gives some indication of 
the overall versatility of POST. 


TABLE 5.- TYPICAL APPLICATIONS OF POST 


Type of 
Mission 


Optimization 

Variables 

Typical Constraints 

CPU Time Required 
to Solve Problems, 
rain 

Type of Vehicle 

Equality 

Inequality 

Asctnc to 

Near-Earth 

Orbit 

Titan II1C A DAE , Space 
Shuttle, Single Stage 
to Orbit (VTO and HTO) 

Payload, Weight at 
Burnout Fuel, Bumtime, 
Ideal Velocity, 

Initial Weight 

Radius 

Flight Path Angle 
Velocity 

Dynamic Pressure 
Accelerations 

2 - 20 

Ascsnt to 
Synchronous 
Equatorial 
Orbit 

Titan I1IC, Shuttle/ 
Tug 

Payload 

Apogee 

Perigee 

Inclination 

Dynanlc Pressure 
Angle of Attack 
Pitch Races 

3-50 

Ascsnt 

Abort 

Space Shuttle 

Abort Interval 

Landing Site 
Longitude and 
Latitude 

Acceleration 
Dynamic Pressure 

2 * 5 

I CBM 

Ballistic 
Miss 11a 

Titan II, 
Klnutaaan I 4 11, 
Safeguard 

Payload 
Mlac Distance 

Latitude 
Longitude 
Cross range 
Down range 

Flight Path 
Angle at Entry 
Acceleration 
during Entry 

2 - 20 

Raantry 

Space Shuttle, X-24C , 
Single Stage to Orbit 

Heat Rate 

Total Heat 
Crossrange 

Latitude 
Longitude 
Crossrange 
Down range 

Heat Rate 
Acceleration 

3-15 

Orbital 

Maneuvers 

Trans tags. 

Space Tug, XUS, 
Solar Electrical 
Propulsion 

Payload 

Fuel 

Radius 

Velocity 

Flight Path Angle 
Argument of 
Perigee Period 
Apogee, Perigee 

Attitude Angles 
Perigee Altitude 


Aircraft 

Performance 

X-24B and C, Space 
Shuttle Subscale, 
Subsonic 

Jet Cruise, Hypersonic 
Bombers and 
Interceptors 

Mach Number 
Cruise Time 
Payload 

Down range 
Crossrange 

Velocity and 
Mach Altitude 

Dynamic Pressure 

Dynamic Pressure 
at Max Altitude 

0.1 - 5 
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The generality of POST that makes it ideal for detailed vehicle per- 
formance work also gives it a unique capability for estimating performance for 
advanced vehicle concepts. A few applications of POST to advanced vehicle/ 
mission concepts are (1) single-stage to orbit trajectory optimization, (2) 
hypersonic cruise vehicle optimization, (3) parachute simulations for solid 
rocket booster recovery concepts, (4) guidance algorithm development for ma- 
neuverable reentry vehicles, (5) ascent optimization for a Mars sample return 
mission, and (6) simulation of a Mars entry using advanced video guidance con- 
cepts . 


PROGRAM STRUCTURE AND FORMULATION 


Executive Design 

POST executive logic was designed to provide (1) readily learned input and 
output procedures, (2) flexible trajectory simulation, (3) detailed vehicle 
modeling, and (4) generalized targeting and optimization. These goals are met 
by the POST executive structure, which is presented in Figure 8. As indicated, 
executive routines are used throughout. These routines control the program 
execution flow by calling subroutines containing the actual mathematical com- 
putations. This modular structure allows the program to be modified quickly 
and easily. 

Both 3D POST and 6D POST are structured in three overlay levels, (0, 0), 
(1, 0), and (2, 0), respectively. The first overlay (0, 0) is the master 
executive overlay, which controls the overall program. This overlay controls 
the read-in of input data and determines which trajectory computations are to 
be performed. All utility routines are contained overlay in (0, 0). 

Overlay (0, 0) first calls overlay (1, 0), which reads the namelist input 
data from cards and stores the processed data on discs for later use. 

Overlay (2, 0) is called by (0, 0) after (1, 0) has completed the input 
processing tasks. The first decision in overlay (2, 0) concerns the type of 
simulation; i.e., single trajectory or search/optimization mode. If a single 
trajectory is to be run, the program calls overlays (2, 1) (2, 2), and (2, 3) 
sequentially, then returns to the master overlay (0, 0). If the search/optimi- 
zation mode is to be used, the program control is turned over to subroutine 
MINMYS, which calls overlays (2, 1), (2, 2), (2, 3), (2, 5), and (2, 6) as re- 
quired to perform the search/optimization function. When convergence has been 
achieved or the maximum number of iterations has been exceeded, control reverts 
back to the master overlay (0, 0) for the next problem. 
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An outline of the approximate calling sequence for each routine is included 
in Figure 8, and shows which subroutines are called by a given routine, thereby 
allowing the detailed logic flow to be followed easily. The overall program 
logic described by the overlays is as follows: 

1) Overlay (2, 1) reads the previously processed input data from tape, 

locates the data for the current phase (event) , and initializes the 

program values based on this input; 

2) Overlay (2, 2) initializes the equations of motion for the current 
phase; 

3) Overlay (2, 3) integrates the equations of motion from time t^ to a 
specified stopping condition for the current phase; 

4) Overlay (2, 5) calculates the control corrections based on the search/ 
optimization algorithm being used, limits the control parameters that 
violate the control parameter constraints, and tests for convergence; 

5) Overlay (2, 6) prints out an iteration summary at the end of each iter- 
ation. It also performs any other information output tasks required by 

search/optimization algorithm, such as printing trail step summaries. 

The program dictionary (subroutine DICT) performs a one-to-one mapping of 
variables in common and the Hollerith names by which the user can select the 
variables for a variety of uses, including output, stopping conditions, control 
variables, and targeting variables. 

POST uses a generalized table storing and lookup procedure whereby the size 
of tables is limited only by the total data storage allocation of 1500 cells. 

Each table has its own multiplier. This is accomplished by dimensioning the 
table by (2) . The first location contains the address of the table and the 
second location contains the table function multiplier. The generalized table 
lookup (GENTAB) is set up to handle all allowable types of tables, namely, 
constant-value, monovariant, bivariant, and trivariant. 

The interaction between the targeting/optimization logic and the trajectory 
simulation was designed so that the trajectory calculations represent an ex- 
ternal evaluation of the objective function and the constraints. In this manner, 
any change to the trajectory simulation capability of the program automatically 
is available to the targeting and optimization logic. 

Coordinate Systems 

POST uses numerous coordinate systems to provide the necessary reference 
systems for calculating required and optional data. The key computational coor- 
dinate systems are illustrated in Figure 9 and 10, and the definitions of all 
coordinate systems used are summarized in Table 6. The equations of motion are 
integrated in the Earth-centered Inertial (ECI) frame. Thus, the external 
thrust and aerodynamic forces, which are computed in the body coordinate system. 
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TABLE 6.- POST COORDINATE SYSTEMS 


Coordinate Systsm 


Axes 


Earth-Cantsred Inartial (ECI) /x^. 


V 


Earth-Centered Rotating (ECR) 


(V V 


Definition 

Earth-centered Cartesian system with coincident with the North 
Pole, Xj coincident with the Greenwich Meridian at time zero and 
and in the equatorial plane, and y^ completing a right-hand 

system. The translational equations of motion are solved in this 
system. 

z R 1 Similar to the ECI system except that it rotates with the Earth 

' so that x^ is always coincident with the Greenwich Meridian. 


Earth Position Coordinates 


Geographic (G) Axes 


Inertial Launch (L) Axes 


Body Reference (BR) Axes 


Orbital Elements 


Vernal Equinox (VE) Axes 


Body (B) Axes 


(V e> h ) 

(v V 2 g) 

(*L * y L * *t) 


These are the familiar latitude, longitude, and altitude desig- 
nators. Latitude is positive in the Northern Hemisphere. Longi- 
tude is measured positive EaBt of Greenwich. Altitude Is measured 
positive above the surface of the planet. 

This system is located at the surface of the planet at the 
vehicle’s current geocentric latitude and longitude. The x^ axis 

is in the local horizontal plane and points North, the y^ axis is 
in the local horizontal plane and points East, and z^ completes a 
right-hand system. 

An Inertial Cartesian system that la used as an inertial reference 
system from which the Inertial attitude angles of the vehicle are 
measured. This coordinate system is automatically located at the 
geodetic $ and inertial longitude of the vehicle at the start of 
the simulation. 


( X BR* y BR* Sr) 


(V V *< “• 9 - n ) 


(\T y VE' z ve) 


Right-hand Cartesian system aligned with the body axes as follows. 
The axis is directed along the negative Xg axis, the y^ R axis 

is directed along the positive axis, and the Zg^ is directed 
along the negative Zg axis. 

A nonrectangular coordinate system used in describing orbital 
motion. The orbital elements are apogee altitude, perigee alti- 
tude, inclination, longitude of the ascending node, true anomaly, 
and argument of perigee. 

A 1950 mean equator and equinox Earth centered inertial system. 

The Xyg axis is in the equatorial plane and is directed forward 

of the vernal equinox of 1950, the z^ axis is directed along the 
north pole, and y^^ completes the right hand system. 

The body axes form a right-hand Cartesian system aligned with the 
axes of the vehicle and centered at the vehicle's center of 
gravity. The Xg axis is directed forward along the longitudinal 

axis of the vehicle, y^ points right (out the right wing) , and Zg 
points downward, completing a right-hand system. 


must be transformed to the ECI system. This is performed by the inverse of the 
transformation matrix, [IB]. The [IB] matrix is functionally dependent on the 
attitude of the vehicle, and is calculated based on the equations describing the 
attitude steering option selected by the user. POST contains four standard 
attitude reference systems as described in Figure 11. Any given trajectory 
problem may require use of one or more of these systems. For example, simula- 
tion of a complete ICBM trajectory typically involves the use of inertial Euler 
angles during first and second stage boost, relative Euler angles during third 
stage flight, inertial aerodynamic angles during the postboost maneuvers, and 
finally aerodynamic angles during reentry. The availability of all these op- 
tions, while confusing to new users, is a valuable aid to the experienced tra- 
jectory analyst, and enables him to steer each phase of the trajectory in the 
most appropriate manner. 
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The matrix transformations between each pair of the coordinate systems 
are presented in Volume I. The basic relationships between the principal com- 
putational systems are depicted in Figure 12. The inverse transformations be- 
tween these coordinate systems can be easily computed by merely transposing 
the matrix elements because of the orthonormality of these matrices. 


Planet Model 


The planet model consists of three basic categories of equations and input 
data: (1) oblate planet geometry and constants, (2) a gravitational model and 

its constants, and (3) an atmosphere model that includes winds. In each of 
these categories, the user can select prestored options to minimize the amount 
of required input data. On the other hand, if the desired option is not avail- 
able, then the user can define his own model via input constants. The only 
inherent program limitations are defined by the equations representing the vari- 
ous models. 


The oblate spheroid model is defined in terms of the equatorial radius R^,, 

the polar radius R , the rotation rate fl , the gravitational constant y, and 
p P 

the second, third, and fourth gravitational harmonics, J 2 , d 3 , and J 4 , respec- 
tively. The 1960 Fisher Earth models are preloaded in POST. The geometry of 
this spheroid is Illustrated in Figure 13. 



Figure 12.- Transformations 

between Coordinate 
Systems 


North Pole 



Vehicle 


Figure 13.- Oblate Planet Geometry 


The gravitational model includes optionally second, third, and fourth 
harmonic terms. This model is adequate for near— Earth ascent, on— orbit, and 
entry performance work. For extremely high altitude satellite maneuvers or 
ephemeris prediction, a more detailed model that includes lunar and solar per- 
turbations and higher order harmonics is generally required. These terms are 
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not included in the standard POST program but are available in a special 
orbital version. 


The gravitational accelerations calculated are given by the equations 


G XI = ?r Hz, r) 


where 


3 XI “ “ 7T p < z > r > 


ZI 


H 

r3 


(l + JR 2 (3 - 5Z 2 )) 

2 + H — ( 

' / 

r \ 


+ DR 4 


- 10Z 2 + 9Z 4 J z 


P(z, r) 


1 + JR 2 (1 - 5Z 2 ) + H ~ (3 - 7Z 2 ) z + DR 4 ^9Z 4 - 6Z 2 + 


and 


(i) 


x * Xj* y ■ y r . z - 2 j, r = r v 

R " Vh* Z = Z l/ r i* J = | J 2. H “ f J3. D = - f- J 4 

POST has the optional capability of three atmospheric models — the general 
table lookup, the 1962 U.S. Standard atmosphere and the 1963 Patrick AFB atmos- 
phere using polynomials. The general table lookup model gives the user the flexi- 
bility of inputting his own atmospheric model if none of the preloaded models is 
adequate. This is particularly useful in performing trajectory analysis for 
planets other than Earth. 

The table lookup atmosphere model can be defined entirely by using tables 
that show pressure, temperature, speed of sound, and density as functions of 
altitude. The speed of sound and density tables can be omitted if desired; in 
this case, the speed of sound and density are computed as 

c s --fir* 


P 


= K 2 


P 

T 


( 2 ) 


yR* 0 

where ® r: » K 2 ® y is the ratio of specific heats, M is the molecular 

0 u 
weight, and R* is the universal gas constant. The equations and constants de- 
fining the 1962 U.S. Standard and the 1963 Patrick AFB atmosphere models are 
given in reference 4. 
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The atmospheric wind velocity components can be input as tables using 
either meteorological or vector notation. If these tables, which are normally 
functions of oblate altitude, are not input, then the atmosphere is assumed to 
rotate uniformly with the planet. 


The wind velocity components can be input directly in the geographic frame 
by defining u w> v w , and w w> or by defining the wind speed (v w j , the wind azimuth 

and the wind azimuth bias . The resulting wind velocity components 

in the G-frame are: 


V 


V w (h) cos (A zw (h) + A zwe ) 


V h > sin ( A zw (h) + ^wb) 

w w (h) 


(3) 


It is clear from the above equation that to input vector wind data A zwg must be 

Input as zero, whereas for meteorologic data the preloaded value of 180 deg 
should be used. 


The wind velocity in the ECI frame is then given by 


^WI 


[IG] 


-1 




Thus, the atmospheric relative velocity vector in the ECI frame is 




( 4 ) 


Vehicle Model 


The vehicle model consists of five general categories of equations and 
data: (1) mass properties, (2) propulsion, (3) aerodynamics, (4) aeroheating, 

and (5) guidance, navigation, and control (GN&C) . The mass properties model 
includes the calculation of vehicle and propellant weights, moments and products 
of inertia, and the location of the vehicle center of gravity. The propulsion 
model computes engine thrust and flowrate for as many as 15 separate engines. 
Standard equations for modeling rocket, turbojet, and ramjet engines are avail- 
able in the program. The aerodynamic model includes all standard ways used to 
describe aerodynamic force and moment coefficients for both 3D0F and 6D0F work. 
The aeroheating model includes all popular heatrate calculation methods, such 
as the standard Chapman's equation for laminar flow and a nonstandard maximum 
centerline technique developed for Space Shuttle. The GN&C model includes all 
equations used for 3D0F openloop steering, as well as general modules for 
specific 6 DOF simulation of an actual flight control system. These modules 
include: (1) a sensor module, (2) a navigation module, (3) a guidance module, 

(4) an autopilot module, (5) a controls module, and (6) an airframe module. 
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Several specific models are available in each of these general modules; however, 
in most cases these modules must be either modified or replaced when a new 
system is to be simulated. Procedures for making these modifications or addi- 
tions are described in reference 5. 

Mass properties .- There are numerous options available for specifying the 
initial gross weight of the vehicle, calculating the time rate of change of the 
vehicle's gross weight, and computing the amount of weight to be jettisoned at 
specified events. The details of these options are presented in Volume II, and 
only the general principles are given in this document. 

The basic approach employed to compute the time history of the vehicle's 
gross weight is to specify only the gross weight at the first event and then 
subtract weight losses due to propellant consumption, ablation, and staging. 
Using this approach, the gross weight of the vehicle at the beginning of the 
simulation is given by 


VT 


+ 

ST G 


W. 


PLD 


( 5 ) 


where W gTG is the gross weight without payload and W pLD is the payload weight. 
The weight of the vehicle during any particular phase is given by 

W G (t) = W G “/*G dt <*) 

where W G is the gross weight on the positive side of the event defining the 
beginning of the phase, and W G is the total derivative of gross weight result- 
ing from engine flowrates and/or thermal protection system ablation. For events 
other than the first, the change in gross weight across the events is computed 
as 


wt = w“ - w, - w 

G G jett PR 


( 7 ) 


where is the jettison weight, and W pR is the weight of propellant remain- 

ing in the previous stage. The propellant remaining can be computed for all 
engines or for a single selected engine, and is given by 


W PR = W p - W PC 


( 8 ) 


where W p is the initial weight of propellant and W pc is the amount of propellant 


consumed. This latter term is computed as 

W PC 


w. 


PC 


dt 


( 9 ) 
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where W is total flowrate for the selected set of engines, and W pc is the in- 
itial value used to carry this parameter over a selected set of phases. When 
is set to zero at the beginning of a phase, then W p ^ represents only the 
amount of propellant used during that phase. The jettison weight, W^^^, can 

be computed as an input constant or determined from an input mass fraction 
table. When mass-fraction is used to determine jettison weight, then 


W 


jett 


u p (M 


( 10 ) 


where A is the stage mass fraction computed from a user's input table. The 
composite center of gravity and the inertia matrix are input in the vehicle 
reference system as defined in Figure 14. In 6D POST, the moments and pro- 
duces of inertia are defined as the integrals 
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/y 2 + z 2 dv 
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zz 


/x 2 + y 2 dv 
and the inertia matrix is given by 
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( 12 ) 


Generally, the center of gravity coordinates and the elements of the inertia 
matrix are input as tables with gross weight as the independent variable. 

Propulsion.- POST can simulate both rocket and jet engines. As many as 
15 engines can be used in either mode in any simulation phase. The equations 
used to calculate net thrust and flowrate per engine are summarized in Table 7. 
The thrust equation, computed in the Body frame, is given by 

p - T u (13) 

— TB ± i^L 

where T^ is the net thrust of the i-th engine, and u^ is a unit vector along 

the thrust direction. In the general case, the direction of the thrust is 
computed as 
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cos 6e cos 6e 

p y 


sin fie 

y 


- cos fie sin fie 
y P 

where fie p and 6e y are the pitch and yaw engine gimbal angles as defined in 

Figure 15. In most 3D0F simulation work, u is simply a unit vector in the 
direction of the x-body axis. 



1) Moments and Products of Inertia in Body Frame (Tables) 

2) Engine Gimbal Points in Body Reference Frame 

3) Aerodynamic Reference Point in Body Reference Frame (Tables) 

4) Center of Mass Location in Body Reference Frame (Tables) 

Figure 14,- Mass Properties Input 
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TABLE 7.- BASIC THRUST AND FLOWRATE EQUATIONS 




Figure 15.- Engine Gimbal 
Angles 

Aerodynamics . - The aerodynamic force coefficients can be expressed in 
terms of the axial force, normal force, and side force, C^, C^, and C^, re- 
spectively. Here and produce forces that act in the -x fi and -z^ direc- 
tions, and Gy produces a force acting along y fi . The aerodynamic force coeffi- 
cients can also be expressed in terms of the .lift, drag, and side-force 
coefficients C , C , and C (Figure 16), where C and C are directed normal 

Li U X Li D 
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to, add along the velocity projection, in the Xg-Zg plane. Note that C y pro- 
duces a side-force, F , acting in the direction of y 


s 



Lift and drag force coefficients are transformed internally to axial and 
normal force coefficients as follows: 


N 


cos a -sin a 

sin a cos a 


(15) 


where a is the angle-of-attack. 


The aerodynamic moments coefficients are defined as: - roll, C 

£ m 

pitch, and - yaw. The pitch and yaw moment coefficients are used in the 3D 

POST static trim option, and the roll coefficient is added only in the 6 DOF 
version. 


All aerodynamic coefficients can be input as constant, monovariant, bi- 
variant, or trivariant tables. In general, there are four tables allocated to 
each coefficient in 3D POST, and eight tables per coefficient in 6D POST. Most 
of these tables can have arbitrary mnemonic multipliers. The mnemonic multi- 
plier capability enables either the coefficient or its slope to be input di- 
rectly into the program. The mnemonic multipliers can be input as the name of 
any computed variable in the output variable list. The coefficient for a given 
table is then the value of the table lookup multiplied by the value of the 
variable defined as the mnemonic multiplier. The values of the aerodynamic 
force coefficients are computed by summing individual contributions as defined 


28 



by user input. As an example, the axial force coefficient is computed from the 
general equation 


C A ■ \ + C A + 


r C p jp + C y «y — 3D POST 
A 5 A 6 


C. a 6 a + C A e 6 a 4* C A r 5r + ^ C f 6 f — 

> A s A « A s w 6 


6 D POST 


(16) 


where each of the coefficients C through C f 3 can be defined as functions 

A 0 A 6 

of as many as three variables, and the mnemonic multipliers through 6 f 3 
can be defined by Hollerith input. Typical examples of are: 


C 


A 


- C A (M) + C a (M) cx 

X° x 

monovariant tables 


mnemonic multiplier 


or 


C. = 


C A (a,M) 


\ 


X 


no mnemonic multiplier 
bivariant table 


In 3D POST, 5p and 6 y are generalized pitch and yaw control surface deflection 
angles. Similarly, in 6 D POST 6 f , i * 1,2,3, are general deflection angles, 

and 6 a, 6 e, and 6 r are the aileron, elevator, and the rudder deflection angles, 
respectively. The detailed equations for all aerodynamic coefficients are 
given in the formulation manuals. 



(17) 


where the dynamic pressure q is given by 



and S is the aerodynamic reference area. 
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Aeroheating .- POST provides for a wide variety of aeroheating calculations. 
Some of these options are specific in nature and apply only to particular ve- 
hicles, whereas others are quite general. The general option is based on tri- 
variant table interpolation of a user-defined heat rate table and provides com- 
plete flexibility with regard to vehicle shape and heat transfer methods. An 
approximate maximum centerline heating rate option is also provided. This 
approximation is based on an analysis of heat rate data calculated by the 
MINIVER aerodynamic heating computer program. A least squares curvefit was 
used to obtain fifth order equations that describes the curves for reference 
heat rate and for the altitude-velocity and angle of attack corrections. The 
equations in this option were written for centerline locations aft of the nose 
of the vehicle. These equations are presented in the Formulation Manual. The 
other heat rate models are based on Chapman's equation 



where the scale factor, K, can be computed from a table look-up. 

In addition to the basic heat rate calculation, POST contains several aero- 
dynamic heating indicators that provide useful information associated with the 
heating environment. One such indicator for zero total angle of attack is 

t 

«I ’/ " V A dt - 

o 

which can be modified for various nonzero angle of attack situations. 

These heat rate indicators can be used in conjunction with a ten panel 
vehicle heating model to incorporate the heating calculation in the vehicle 
weight calculations. 

In this simplified model, the total heat for each panel is assumed to be 
a constant ratio of the total heat calculated by the selected option at a given 
location on the vehicle. The total thermal protection system (TPS) weight is 
then computed as the sum of the individual panel weights, 

i 

10 

W TPS = W uA A i 
i=l 1 


where W 


uA. 


is the weight per unit area and 


is the area of the i-th panel. 
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Guidance and Control.- POST can simulate both open loop and closed loop 
guidance. In 3D POST, the vehicle's actual attitude is computed directly from 
an attitude equation (polynomial, table, etc.) or from a user-programmed guid- 
ance equation. This approach simulates a perfect autopilot where the actual 
and the commanded attitudes are identical. In 6D POST, the complete flight 
control system can be simulated in substantial depth. However, this generally 
requires that the user code the detailed models associated with the particular 
sensor, navigation, guidance, autopilot, and control systems being simulated. 


The sensor module computes information that describes the behavior of the 
sensing elements of the vehicle's navigation system. Thus, the primary func- 
tional responsibility of the module ia that of simulating hardware character- 
istics of sensors. For example, the behavior of an inertial measurement unit 
(IMU) can be described by a mathematical model of the platform and the accel- 
erometers. Frequently this module is used for error analysis purposes. 

Sensor models called by this module are necessarily vehicle and subsystem 
dependent. As a consequence, the sensor model must be designed and imple- 
mented for each particular application. 

There are many applications of the program that do not require a specific 
simulation of the sensors. Therefore, for convenience, a "perfect" sensor 
model is coded into this routine. This "perfect" sensor model sets the sensed 
program variables equal to their actual values as calculated in the simulation 
models . 

The function of the navigation model is to estimate the state (position, 
velocity, etc.) of the vehicle based on the sensor outputs. Clearly, this 
module is also vehicle and subsystem dependent and must be designed and imple- 
mented for each specific application. The 6D POST contains no built-in navi- 
gation models. As a consequence, the estimated state is set equal to the actual 
state. This is also equivalent to simulation of perfect navigation. 

The guidance module takes the output of the navigation model and computes 
a guidance command. Typically, the guidance command represents a desired change 
in the current attitude of the vehicle. This command is computed on the basis 
of meeting some specified trajectory condition, such as an inject condition or 
a landing condition. The autopilot is designed to remove the errors between 
the commanded values of the guidance variables and their actual (or sensed) 
values. This is accomplished by deflecting engines, control surfaces, and/or 
firing RCS jets. 

The current version of 6D POST contains three preloaded guidance options: 
(1) an open-loop profile steering; (2) a closed loop v-h profile ascent algo- 
rithm; and (3) the constant drag Space Shuttle reentry scheme. If these methods 
are inadequate, the user may implement his own guidance algorithm into this 
module . 
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The function of the autopilot module is to generate a command that, when 
implemented through the deflection equations contained in the controls module, 
causes the vehicle to respond as prescribed by the guidance module. This func- 
tional responsibility is depicted in Figure 17. 



L — .. . from Simulator 

Nomenclature : 

6 - Guidance command 

— c 

A a - Actual or sensed values of guidance variables 

e - Generalize error signal 

69. ~ Pitch, yaw, and roll autopilot command 

A - Deflection (engine or aerodynamic surfaces) angles 

Figure 17.- Guidance and Control Block Diagram 


The autopilot module calculates only autopilot commands based on the input 
guidance commands, and does not calculate engine or control surface deflections. 
The engine and control surface deflections are computed in the controls module 
as a linear function of the autopilot commands. The autopikot commands 60, 6$, 
Sip represent changes in vehicle attitude. The mixing equations determine the 
engine and control surface deflections that create the control forces and 
moment s . 

Currently, there are two Space Shuttle autopilot models available in 6D 
POST. One autopilot is for ascent and the other for reentry. The ascent auto- 
pffoi fs somewhat standard and could be used on most ascent problems with little 
or no modification. The basic inputs to this model are: attitude commands 

from the guidance, inertial attitude angles, body rotational rates, transla- 
tional accelerations, and preloaded engine deflection commands. The outputs 
are pitch, yaw, and roll autopilot commands, which are sent to the controls 
module to determine the engine deflection angles. The reentry autopilot is 
Space Shuttle-oriented and is probably not applicable to other vehicle 
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configurations. This model is intended to provide attitude control for Space 
Shuttle beginning at an altitude of approximately 400 000 ft and ending in the 
high subsonic flight regime. The control logic makes use of both aerodynamic 
control surfaces torques and reaction control jets. 

The controls model converts pitch, yaw, and roll autopilot commands to 
aerodynamic control surface deflection angles and/or engine gimbal angles. The 
conversion of the autopilot commands into deflection angles is implemented 
through the matrix mixing logic given by the equation 

6*6 + [M] 66 (18) 

— o 

where 6_ denotes a general deflection angle with a null position of 6^, [M] the 

mixing gains** and 60_ the autopilot commands. The gains contained in the mixing 
matrix, [M], and the null deflections, 6^, are specified by user input. 

In 3D POST, there are five basic types of openloop guidance options for 
controlling the attitude of the vehicle during 3 DOF trajectory simulation. 

These options are as follows: 

1) Body rates; 

2) Aerodynamic angles; 

3) Inertial aerodynamic angles; 

4) Relative Euler angles; 

5) Inertial Euler angles. 

The body rate option is generally used to simulate strapdown systems with the 
body rates being computed from user-specified polynomials. The attitude of the 
vehicle is then determined by integration of the quaternary equation 

e =• \ [E] u. (19) 

When using this option the user must specify (1) the initial attitude of the 
vehicle, and (2) the coefficients and the Hollerith names of the arguments of 
the polynomials used to compute the body rates. 

Atmospheric and inertial velocity relative aerodynamic angles are gener- 
ally used to simulate reentry and orbital maneuvers, respectively. Similarly, 
relative and inertial Euler angles are typically used to simulate vehicles that 
employ local horizontal or inertial reference systems. In all of these appli- 
cations, the attitude angles can be computed based on five basic techniques: 
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1) Polynomial steering: Under this option the steering angles are com- 

puted from a cubic polynomial 


3 

0 ( y ) = c^y 1 (20) 

i=0 

where the user selects the coefficient c and the independent vari- 
able y. The highest-order coefficient that is input determines the 
degree of the polynomial. The argument can be selected as any in- 
ternally computed variable; e.g., time, velocity, altitude, etc. The 
constant term of the polynomial, c q , can be either input at the be- 
ginning of each phase or carried across as the value of the angle at 
the end of the previous phase. The polynomial coefficients are gen- 
erally used as the independent variables for targeting/optimization. 

2) Table steering: Under this option the steering angles are computed 

via table interpolation, which is denoted by 

6(X) - e m T n [f( Z )]. (21) 

The user initially inputs the table multiplier 0 , the order of 

m 

interpolation n, and the table data [^, f(^)]. The table multiplier 
or the dependent table values can be used as independent variables 
for targeting or optimization. The order of interpolation can either 
be linear or cubic. The tables* can be monovariant, bivariant, or 
trivariant functions of the table arguments. 

3) Piecewise linear steering: Under this option the steering angles are 

computed from a general piecewise linear function of the form 


B(y) = ci + 


C 2 ~ Cl 

Y2 - yi 


(y - yi) 


( 22 ) 


where C] is equal to 0 at the beginning of the current phase, C 2 is 
the desired value of 0 at a designated later event, is equal to 

the value of the designated event criterion at the beginning of the 
current phase, y 2 is the desired value of the designated event, and 
y is the current value of the designated event criterion. 


This option is similar to the polynomial option except that the values 
of 0 are specified directly rather than as 0 q, 0, 0 and Clearly, 

0 is linear in time if y = t; otherwise 0 is only linear in y. When 
the desired values of the steering angles are used as independent 
variables, the problem of cascaded steering effects is avoided and 
the target ing/op timization algorithm generally converges faster. This 
option also automatically computes the steering angle rates required 
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to change the attitude to the desired value at the designated event, 
which reduces the problems related to guessing accurate initial pitch 
rates . 

4) Linear feedback steering: Under this option the steering angles are 

computed from the linear error-error rate feedback control law 

6 - c i + K D ( F a - F d) + k r(l - *d) (23) 

where c^ is a nominal steering angle profile, is the displacement 

error gain, K is the rate gain, F - F is the error in the steering 
R 3 Q 

function error. 

This option is, of course, the classic path control law, and enables 
the user to steer to a wide variety of trajectory profiles, such as 
velocity versus altitude profile, acceleration versus altitude pro- 
file, etc. This option is particularly useful for reentry trajectory 
shaping. 

5) Generalized acceleration steering: Under this option the steering 

angles are computed by solving a set of user-specified equations. The 
dependent and independent variables in these equations must be selected 
from the dictionary of variables already computed in POST. The only 
restriction is that these equations must be explicitly a function of 
some derivative -compound in the inner loop of the program. 

In more precise terms, the steering variables are determined from the 
iterative solution of the problem: 

For each instant of time, determine the values of the steering vari- 
ables, 6, that satisfy the steering equations, 

e(0) = y_(0_) - Zrf = 0 (24) 

where is a n— component vector of dependent variables , y j is the de- 
sired value of these variables and e the error in dependent variables. 

A typical application of this option is control normal acceleration 
to one-g and axial acceleration to three-g by calculating the angle 
of attack and throttle setting that solves the equations 

A^Cct, n) - 96.6 = 0 

A z B (ct » n) - 32.2 = 9 (25) 
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Trajectory Simulation 


The trajectory simulation features consist of the integration of both the 
translational and the rotational equations of motion, the event interrupt and 
sequencing logic, and the numerical (in some cases analytical integration) 
methods. The flexibility that the user has in selecting these three features 
is the key to the utility of the POST programs. 

Translational equations .- The translational equations of motion are solved 
in the Earth-centered inertial coordinate system (ECI) . These equations are 


j-I + W' 1 [*TB + 4 b] + £i (26) 

where is the thrust acceleration in the body frame, A ^ is the aerodynamic 

acceleration in the body frame, and is the gravitational acceleration in 

the ECI frame. The external accelerations are first computed in the body frame 
and then transformed to the ECI frame. The external forces that cause these 
accelerations are illustrated in Figure 18. In 6D POST, the net translational 
force due to the reaction control system are also included in the total non- 
gravitational force acting on the vehicle. 



mVj- [IB]" 1 ^ + FJ + mGj; 


Figure 18.- Force System in Pitch Plane 
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There are a number of options for initializing the translational equations 
of motion. These options are summarized in Table 8. The most frequently used 
option is the planet relative spherical coordinates (<t>, 0, h) , and the local 
horizontal velocity components /V R , y R > A zr) * 


TABLE 8.- TRAJECTORY INITIALIZATION OPTIONS 


Position 

Velocity 

Attitude 

Angular Velocity 

Inertial Rectangular 

Inertial Rectangular 

Inertial Euler Angles 

Body Rates u yt 

(*i* y i* z i) 

( V *r V yl * v zi) 

(V V 9 i) 

Inertial Euler Rates 

Earth - Relative Polar 

Inertial Local Horizontal 

Relative Euler Angles 

/i-. a . 

(r, ♦, 8) 

( V I* Y I* A zi) 

r*’ v * R ) 

( ¥ i* i’ i) 

Orbital Element 

Earth-Relative Local 

Aerodynamic Angles 

Relative Euler Rates 

(a, e, i, n, u)» 6) 

Horizontal (V R , y r , A zr ) 

Atmospheric Relative Local 
Horizontal | V^, y a * 

Orbital Elements 
(a, e , if ui , 6) 

(ot, 3, cr) 

(V V M 


Rotational equations .- In 6D POST, the rotational equations of motion are 
solved in the body-centered coordinate system. These equations are 

i - 1 [ E] ^ 

“b = W 1 [Mg - [i] “b - % x W “b] 

^B = — TB + — AB U 


where 


N 

eng 

^TB " "2 “TB X A ^BT 
i*l 
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and 


— AB = 


d R C * 


d P C » 


d Y C n 


- — AB X A 4b 



X ref 

- X 

eg 

a ^ab = 

y ref 

- y cg 


_* Z ref 

- z 

eg J 


In the rotational equations, £ is a four-dimensional vector of quaternary 
parameters, [E] is the quaternary matrix, u>g is the inertial angular velocity 

expressed in the body frame; Mg is the total external moment acting in the 

vehicle as a result of the thrust, the RCS, and the aerodynamic forces, and 
II] and [EJ matrices are given by 


[I] » 


I XX “ I XY *XZ 


I XY I YY " I YZ 


X XZ ” I YZ : ZZ 


[E] = 


S 1 e 2 e 3 


®0 ®2 " e 3 


6 0 “ e l e 3 


e e — e 
0 1 e 3 


(28) 


where the inertia matrix, I, is not necessarily assumed to be constant-valued. 
The rotational equations of motion are initialized by defining both the initial 
attitude angles and the initial attitude rates. There are three options for 
tislizat ion as summarized in Table 7. The equations for each of these 
options are presented in reference 6. 
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Events and phases .- The concepts involved with the definition of the tra- 
jectory sequence of events, and the relationships between these events and 
simulation phases are fundamental to the use of POST programs . As a result, 
special emphasis should be placed on clearly understanding these simple yet key 
concepts . 

In POST terminology, an event is defined in terms of three critical ele- 
ments: (1) an event number, (2) event criteria, and (3) an event criteria 
value. These three elements combine to define a condition, which when it 
occurs, caused the trajectory simulation to be interrupted. The ability to 
interrupt the simulation based on any user-defined condition can be used for 
a variety of purposes. For example, an interrupt can be used to change vehicle 
data, such as aerodynamics or propulsion! or to change simulation control data, 
such as integration methods and so on. 

Mathematically, the i-th event criteria defines a scalar-valued continuous 

function, y (t), and the event numbers index and, in most cases, order these 
i 

event criteria. The event criteria value, y ± , is used in conjunction with the 
event criteria, y^t), to define the event interrupt equation 

y^t) - = 0 

The time at which the i-th event occurs is then computed by iteratively solving 

this equation for the value t = t^, where |t^j = y^. This concept is illus- 
trated in Figure 19. 



Figure 19.- Illustration of Time-to-Go Logic 
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Each event (other than the first event) serves to (1) terminate the pre- 
vious phase, and (2) initiate the following phase. This can be explained as 
follows. Event i is said to occur at time t - t^, and the plus (+) and minus 

(-) sides of the i-th event occur at the left and right limit times, t+ and t~, 

respectively. The trajectory data at time t~, and in particular the last 

printblock in each phase, are compyted using the simulation data defined in 

phase i-1. Similarly on the plus side of the event, at t+, the trajectory data 

are recomputed using the new data (if any) defined by input for phase i. This 
subtle concept is very important in the simulations where mass is discontinuous 
due to staging, or accelerations are discontinuous because of changes in pro- 
pulsion and configuration characteristics. 

Four types of events have been defined to provide flexibility in setting 
up a given problem: 

1) Primary events: These describe the main sequential events of the tra- 

jectory being simulated. These events must occur, and must occur in 
ascending order according to the event number. Most problems will 
usually be simulated by a series of primary events; 

2) Secondary events: These are events that may or may not occur during 

the specified trajectory segment. Secondary events must occur in 
ascending order during the interval bounded by their primary events. 
The occurrence of a primary event will nullify the secondary events 
associated with the previous primary event if they have not already 
occurred; 

3) Roving primary events: These events can occur any time after the 

occurrence of all primary events with smaller event numbers. They 
can be used to interrupt the trajectory on the specified criterion 
regardless of the state of the trajectory or vehicle. 

4) Repeating roving events: These events are the same as primary roving 

events except the interrupt values are input differently. There are 
two options for criteria value specifications. Option 1: Input the 

initial value, the increment, and the number of times the event is to 
be repeated. Option 2: Input an array of event criteria values. 

The cycling routine monitors as many as ten events at a time, depending on 
the types of events to determine which event is to occur next. 

Multiple events are monitored in the following sequence: 

1) The next primary event is monitored; 

2) As many as nine primary roving events are then monitored, provided 
there are no secondary events. A roving primary event is added to the 
list of those being monitored as soon as the primary event immediately 
preceding that roving event has occurred; 
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3) Next, as many as nine secondary events are monitored, provided there 
are no primary roving events. (Note that caution must be exercised 
when using secondary events because of their nature. Because as many 
as nine secondary events are monitored at a time, any one of those 
nine will occur as soon as its criterion has been met. Because they 
are secondary events, the event that occurs will cancel all secondary 
events with smaller event numbers.); 

4) Finally, a total of nine primary roving and secondary events are moni- 
tored. 

Because the program can only monitor nine events (in addition to the next 
primary event) , the sum of the primary roving events and the secondary events 
in a phase must be less than or equal to nine or a fatal error will result. 

The time-to-go model iterative solves the event criteria equations and 
determines when the events occur during the trajectory simulation. Basically, 
the algorithm checks the values of the criteria being monitored at each inte- 
gration step. If none of the criteria values has bracketed the desired cutoff 
value, then another integration step is taken. If a criterion variable is 
bracketed with the input step size, then a new stepsize is computed equal to 
the predicted time-to-go. 

The predicted time-to-go for each event is computed from the basic secant 
equation 

At* * -Ay^( t) At/^Ay^t + At) - Ay^t)] (30) 


where Ay(t) is the difference between the actual and the desired value of the 
event criterion. If more than one event is bracketed, then the minimum pre- 
dicted time-to-go is used as the integration stepsize. This process is re- 
peated until the criterion value is within the specified tolerance of the 
desired value. Again see Figure 19. 

In POST, the simulation data are input by phase, where each phase is de- 
fined by two events — the event initiating the phase and the event terminating 
the phase. As a result, the user is required to define via events the sequence 
of trajectory phases that specify from beginning to end the problem being simu- 
lated. Physically the data for each phase are located in the data deck be- 
tween the Hollerith specification of the events defining the phase. 

Specific event numbers, which are selected by the user, must satisfy cer- 
tain conditions based on the type of events employed for a given problem. In 
general, the event numbers must be monotonic increasing, but they need not be 
consecutive. 

Integration methods .- There are three general purpose numerical integra- 
tion methods available as automatic program options; (1) standard fourth order 
Runge-Kutta, (2) an eighth order Runge-Kutta, and (3) a variable-step variable- 
order predictor-corrector. The vast majority of users employ the standard 
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Runge-Kutta integrator for atmospheric flight simulation. However, recent 
computational experience with the predictor-corrector method indicate that it 
is substantially more efficient on certain classes of problems, particularly 
those requiring orbital propagation. This method, developed by F. T. Krogh of 
the Jet Propulsion Laboratories, represents the state-of-the-art in numerical 
solution of systems of ordinary differential equations. It includes all of 
the following facilities: 

1) A core integrator for advancing the solution from one uniform step to 
the next consisting of variable order Adams predictor-corrector equa- 
tions requiring the storage of a difference table for only the highest 
ordered derivatives; 

2) A method to start integration with first-order equations and to in- 
crease the order to as high a level as numerical stability permits; 

3) Algorithms for changing the stepsize and updating accordingly the dif- 
ference tables of the highest— order derivatives including appropriate 
smoothing to prevent numerical instability; 

4) Algorithms for deciding when and by how much to change the stepsize 
based on the accuracy requested by the user; 

5) Tests for numerical stability of the predictor-corrector order and 
stepsize tentatively chosen in the context of the current differential 
equation set; 

6) Algorithms for the automatic selection of the core integrator to fit 
the characteristics of the set of differential equations at hand; 

7) An interpolation algorithm for obtaining dependent variable values to 
the user-specified accuracy at values of the independent variable 
different from normal integration steps. 

In addition to these general purpose numerical methods, 3D POST contains 
two special integration methods for propagating orbital trajectories: (1) 

Laplace's Method, and (2) Encke's Method. 

Laplace's method is an iterative technique for propagating the position 
and velocity (in planet-centered coordinates) of a nonthrusting vehicle in 
vacuum during flight over a spherical planet. The technique is based on the 
analytical solution of the two-body equations, and yields the inertial state 
at time t + At as a function of the state at time t. 

The equations used in Laplace's method are: 

JLjCt + At) » fr^t) + gVjCt) 

V : (t + At) - fr^t) + gV^t) (31) 
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The coefficients f, g, f, and g are computed analytically from the change in 
eccentric anomaly during the time period At. The change in eccentric anomaly 
is calculated by solving a special version of Kepler's equation via the Newton- 
Raphson algorithm. 

The Encke method used in 3D POST is modified from the usual Encke tech- 
nique in that it rectifies the reference conic at every integration step and 
does not use the standard Q-series expansion in calculating the gravitational 
increment. 

The Encke method should be used for orbital problems where small perturb- 
ing accelerations, such as the oblateness of the planet, atmospheric drag, or 
solar electric propulsion, must be included in the simulation. Numerical re- 
sults indicate that, for problems involving small perturbations from Keplerian 
motion, Encke’s method is approximately four times faster than Cowell methods, 
which integrate the total acceleration. 

The Encke method determines the total motion by summing the motion due to 
the two-body equations and the motion due to the perturbations to the two-body 
equations. .The position and velocity in the inertial planet -centered system 
at time t + At is given by 

£ (t + At) - r 2 (t + At) + A_r(t + At) 

VjCt + At) * V 2 (t + At) + AV(t + At) (32) 

where _r 2 » V 2 denotes the Keplerian motion computed by Laplace’s equations; 
that is. 


r_ 2 (t + At) - fr 2 (t) + gV 2 (t) 

V 2 (t + t) - fr 2 (t) + gV 2 (t) 

and (Ar, AV) denotes the numerical solution of the differential equations 
Ar - AV 

4V - [[IB]’ 1 [A^ + AyJ + G r ] - % (r 2 + Ar ) 

Ar(t) - AV(t) « 0, 

where £ 2 |£ 2 + ArJ is the two-body acceleration at + A_r. 
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Special Purpose Options 


In 3D POST, there are numerous special purpose options that aid the user 
in modeling special situations that frequently arise in trajectory analysis. 
These options efficiently simulate these special situations by using specific 
equations that are derived on the basis of the assumptions employed. For ex- 
ample, many orbital maneuvers can be adequately modeled as a discontinuity in 
velocity. This special case does not require numerical solutions and is 
modeled using the classic rocket equation. This and other key special purpose 
options are summarized in Table 9. 


Auxiliary Calculations 

In addition to computing the basic variables, POST also computes numerous 
auxiliary variables that are related to (1) conic parameters, (2) range calcu- 
lations, (3) tracking data, (4) analytic impact calculations, (5) velocity 
losses, and (6) velocity margins. Table 10 summarizes the auxiliary variables 
that can be optionally computed. All these auxiliary variables can be used as 
independent variables in the targeting and optimization, as event criteria in 
the simulation, or as output variables of general interest. Special print 
blocks can be requested for conic calculations, tracking station data, and 
velocity loss information. These print blocks are arranged in a convenient 
format and contain all of the auxiliary variables associated with the particu- 
lar type of auxiliary output requested. In some cases, these complete special 
purpose print blocks are not required. If so, the particular output variable 
desired can be added to the normal trajectory print block by following the 
sample procedures as outlined in references 1 and 7. 


Targeting and Optimization 

The targeting and optimization capability in POST provides the trajectory 
analyst with a set of numerical tools that enable him to solve a broad spectrum 
of trajectory design problems. This capability is achieved by first providing 
the user with the option to specify a variety of problem types, such as full 
rang targeting, unconstrained optimization, and constrained optimization in- 
cluding equality and inequality constraints. Second, the trajectory analyst 
is provided the capability to define the details of the problem formulation 
using simple and direct input procedures. These procedures enable him to select 
the optimization variable, the constraints, and independent control variables 
by specifying only their Hollerith names, their event numbers, the constraint 
values and tolerances, and an initial guess for the control variables. Finally, 
the trajectory analyst is provided a library of numerical algorithms that can 
be used to iteratively solve any standard discrete parameter trajectory problem 
formulation. These algorithms include the popular accelerated gradient projec- 
tion method for nonlinear programming formulations, Newton-Raphson for full 
rank targeting formulations, and the conjugate gradient or the Davidon methods 
for unconstrained optimization or constrained optimization employing penalty 
functions. 
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TABLE 9.- SPECIAL PURPOSE OPTIONS 
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TABLE 10.- SUMMARY OF AUXILIARY CALCULATIONS 


Position and Velocity 

Geocentric Radius 

Oblate Atltude 

Inertia Velocity Magnitude 

Planet Relative Velocity 

Atmosphere Relative Velocity 

Inertial Flight Path Angle 

Planet Relative Flight Path Angle 

Atmospheric Relative Flight Path Angle 

Inertial Velocity in Geographic Frame 

Relative Velocity in Geographic Frame 

Atmospheric Velocity in Geographic Frame 

Inertial Azimuth 

Relative Azimuth 

Geocentric Latitude 

Geodetic Latitude 

Inertial Longitude 

Relative Longitude 

Sensed Acceleration in Body Frame 

Sensed Acceleration in ECI Frame 

Conic Variables (Elliptic and Hyperbolic) 

Orbital Energy 

Apogee Radius and Altitude 

Perigee Radius and Altitude 

Semimajor Axis 

Eccentricity 

Inclination 

Period 

Argument of Perigee 
True Anomaly 

Longitude of Ascending Node 
Angular Momentum 
Perigee Latitude 
Perigee Longitude 
Argument of the Vehicle 
Time to Perigee Passage 
Time Since Perigee Passage 
Eccentric Anomaly 
Mean Anomaly 
Perigee Velocity 
Apogee Velocity 

Velocity Required to Circularize 
Circular Velocity 
Longitude of the Vernal Equinox 
Velocity Margin 

Range Calculations 

Dot Product Range 
Crossrange Reference 
Inertial Reference 
Earth Relative Reference 

Auxiliary Attitude Calculations 

Aerodynamic Angles 
Inertial Euler Angles 
Relative Euler Angles 
Inertial Aerodynamic Angles 
Body Rates 

Inertial Euler Angle Rates 
Relative Euler Angle Rates 


Tracking Data 

Slant Range (Vector & Magnitude) 

Elevation Angle 
Azimuth Angle 
Cone Angle 
Clock Angle 
Space Losses 

Analytic Impact Calculations 

Geodetic Latitude of Impact 
Relative Longitude of Impact 
Altitude of Impact 
Position Vector at Impact 
Velocity Vector at Impact 

Velocity Losses 

Drag Loss (Inertial & Relative) 

Thrust Vector Loss (Inertial & Relative) 
Atmospheric Pressure Loss (Inertial & Relative) 
Gravity Loss (Inertial & Relative) 

Ideal Velocity 

Sun-Shadow Calculations 

Position and Velocity in Vernal Equinox System 

Greenwich Hour Angle 

Cone and Clock Angles of Sun Vector 

Shadow Function 

Multiple Vehicles 
Conic Variables 

Position and Velocity Variables 
Accelerations 

Aerodynamic & Aeroheating 

Aerodynamic Coefficients 

Aerodynamic Forces & Moments 

Dynamic Pressure 

Total Angle of Attack 

Dynamic Pressure & Angle of Attack 

Mach Number 

Reynolds Number 

Heat Rate 

Total Heat 

Atmospheric Densigy Temperature , Pressure 
and Speed of Sound 

Propulsion Calculations 

Vacuum Thrust 
Net Thrust 

Total Thrust & Accelerations 
Flowrate (Total) 

Partial Flowrate 
Propellant Consumed 
Propellant Remaining 
Throttle Setting 
Gimbal Angles 
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All the targeting and optimization procedures used in POST are completely 
numerical. This means that all the function evaluations and derivatives used 
by the various algorithms, are computed numerically from only the trajectory 
simulation results. Thus, the user needs only to define a problem formulation 
that is consistent with the vehicle simulation, and the program automatically 
does the rest. This numerical approach, although somewhat slower in execu- 
tion speed than some analytical methods, offers the advantages of "hands off" 
flexibility and a significantly lower program maintenance overhead. 

Problem formulation .- Discrete parameter methods provide the basis of the 
targeting and optimization formulation and solution algorithms employed in 
POST. The use of the discrete parameter approach means that trajectory control 
variables are defined by a finite set of parameters, called independent vari- 
ables, and given values for these variables the trajectory can be uniquely 
computed. A good example of this approach is the representation of the ve- 
hicle’s pitch attitude by a piecewise linear function. In this parameterization, 
the nodes, 0^, and the pitchrates, 0_^, are the independent variables that can 

be varied to steer the vehicle. A near optimal pitch program can be generated 
using these parameters to optimize a selected performance criterion, say pay- 
load, subject to both mission and vehicle constraints. In recent years, this 
discrete parameter approach to trajectory optimization has replaced the more 
complex calculus of variation techniques. This is because the discrete param- 
eter methods are simpler to implement, understand, and use, and, as a result, 
are substantially more reliable. The reason for this is that discrete param- 
eter methods are based on ordinary vector calculus, which is easier for most 
practicing trajectory analysts to understand than the complex procedures re- 
quired in most variational trajectory programs. 

Using discrete parameter concepts, a vast diversity of trajectory optimiza- 
tion and targeting problems can be formulated that have a common mathematical 
structure. First as described above, there is a vector u of control parameters 
which must be selected to define the trajectory 

x(t) * * 2L>H.* t: £o (33) 

where the state variable 3C is typically composed of the components of the ve- 
hicle’s position, velocity, and mass (or weight) and $[•] denotes the numerical 
integration of the equations of motion. Second, for each trajectory problem 
there is a vector of constraints called dependent variables y_(u) , together 
with a vector of constraint limit values, b_. These constraints are calculated 
at a particular user-specified event. Thus, the dependent variables are com- 
puted as functions of the independent variables, ij, from the general relation- 
ship 

Z (u) * Y(u,X(u)) 

where X(u_) is the composite vector of state, conditions at all events where con- 
straints are defined. Thus, 
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X(„) - Y ® 1 [2£-( c 1 ) • ^ V2o] 

where the s£t |t^| * |t:y^(t) K j for fill i| denotes the times et which the 

dependent variable events occur. The number of constraints, n, can be less 
than, equal to, or greater than the number of control parameters, m. In con- 
strained optimization, n < m, excess degrees of freedom exist in the formula- 
tion enabling some performance criteria to be optimized and the constraints 
to be satisfied. In full rank targeting, n - m, no excess degrees of freedom 
exist, and as a result it is possible to satisfy only the constraints. The 
final case, n > m, generally occurs only during the solution process when an 
excess of inequality constraints become active on some intermediate iteration. 
Finally, for each trajectory, there is an objective function F(u) , which is 
calculated in the same manner as the constraints. The object of the problem 
is then to determine the control parameters, u, which are feasible in that all 
constraint parameters fall within their acceptable ranges and optimal in the 
sense that the objective is minimized (or maximized). 

The general discrete parameter trajectory optimization problem is then 
the well known nonlinear programming problem. Symbolically, it is expressed 
as: 

minimize: F(u) 

subject to: y_(u) a b 

where: u is the mxl column matrix of control parameters, 

F is the scalar objective function of the vector of control 
parameters , 

y_ is the nxl column matrix of dependent parameters, 
b is the nxl column matrix of constraint parameter limits, 


a is the nxl column matrix of constraint parameter relations 
_ (each element is the appropriate relation of the triple 
1. or >) . 

Virtually all trajectory targeting and optimization problems can be cast in 
this structure. Specific examples of nonlinear programming formulations for 
ascent, reentry, and orbital maneuvers are presented in the Sample Case Section. 

Algorithm macrologic .- POST uses an accelerated projected gradient algo- 
rithm (PGA) as the basic targeting/optimization technique. PGA is a combination 
of Rosen's projection method for nonlinear programming and Davidon’s variable 
metric method for unconstrained optimization. The program also contains backup 
single-penalty function methods that use steepest descent, conjugate gradients, 
and/or the Davidon method. 
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The accelerated projected gradient algorithm is based on five intuitive 
working principles as follows: 

One -Dimensional Search . - Using cost and constraint function gradient in- 
formation, a direction of search is established. Then a one-dimensional min- 
imization is performed in this direction on an appropriate function. In this 
manner, a difficult multidimensional optimization problem is replaced by a 
sequence of simple one dimensional minimizations. 

Linearized Constraint Correction . - Assume that the current vector of con- 
trol parameters is outside the feasible region. This correction scheme approxi- 
mates the contours of constant constraints as uniformly spaced parallel hyper- 
planes based on their respective gradients and values for the current control 
parameter vector. Using this approximation the smallest correction to the con- 
trol parameters that would satisfy all active constraints is computed; or that 
failing, minimize the sum of the squares of their violations. One -dimensional 
minimization of the sum of squares of the constraint errors is then performed 
along the direction of this correction to obtain the next iterate of control 
parameters . 

Gradient Projection .- Once a feasible vector of control parameters is ob- 
tained, the negative gradient is resolved into two components; one parallel to 
and one normal to the hyperplane tangent to the boundary to the feasible region 
at the current point. A minimization is then performed along the direction of 
the parallel negative gradient component to obtain the next control parameter 
iterate. The function to be minimized in this one dimensional search is the 
fourth basic principle of the algorithm. 

Estimated Net Cost Function .- Because the constraints are nonlinear, the 
tangent plane only coincides with the boundary of the feasible region at the 
point of tangency; hence, a search along the component of the gradient lying in 
the tangent plane will probably terminate at a point external to the feasible 
region. Therefore, the real object of the search should not merely be to find 
the minimum value of the cost function in the search direction. Rather it 
should be to find a unique point along the search ray that yields on correction 
back to the feasible region a new feasible point with the smallest value of the 
cost function. This point is approximately determined by minimizing along the 
parallel component of the gradient the cost function less an estimate of the 
deterioration of the cost function occasioned by correcting back to the feasible 
region. The estimate is based on linearized constraint-correction formulae. 

Gradient Acceleration . - It is widely known that the convergence of uncon- 
strained gradient algorithms can be drastically improved by using gradient in- 
formation from several iterations to estimate the inverse of the Hessian matrix 
of a quadratic form approximating the cost function. In fact for a cost func- 
tion of m control parameters, it can be shown that a Hessian-inverse estimating 
accelerated gradient scheme converges in m iterations and a conventional 
steepest descent algorithm converges only asymptotically. To similarly accel- 
erate the projected gradient algorithm for constrained problems, it is assumed 
that the cost function is a quadratic in m-q variables over the constraint 
boundary. Here q is the number of active constraints defining the boundary. 


49 



Thus convergence should be accelerated to m-q iterations after the set of 
active constraints that determines the feasible region stabilizes. 

For the widely known special cases of the general nonlinear programming 
problem, the accelerated projected gradient algorithm degenerates to the appro- 
priate special purpose state-of-the-art programming procedures. For example, 
if no constraints are present, the algorithm degenerates to the Davidon-de- 
flected gradient procedure. This procedure has long been considered the method 
of choice for solving unconstrained parameter minimization problems. At the 
other extreme, if the problem has more active constraints than controls, the 
algorithm reduces to the Gauss 1 least squares procedure for minimizing con- 
straint violation. This technique is generally conceded to be the best avail- 
able for solving over-determined systems of equations. Similarly, if the 
number of constraints is precisely equal to the number of controls, the algo- 
rithm becomes the widely known Newton-Raphson procedure for solving systems of 
nonlinear equations. This scheme is certainly the simplest of the efficient 
methods for solving fully determined systems of equations. 

A summary of macrologic for the gradient projection algorithm is pre- 
sented in Figure 20. As indicated, if the initial guess violates any con- 
straints, then the algorithm attempts to satisfy all the constraints by taking 
constraint restoration steps. Once a feasible control parameter vector has been 
found, the algorithm generates a sequence of iteration pairs. Each pair con- 
sists of an optimization step followed by a constraint step. If the user's 
initial control-parameter estimate is not feasible, however, a steadily im- 
proving sequence of constraint-correction steps is undertaken until a feasible 
solution is found. Furthermore, the subsequent optimization step is omitted 
after any constraint-correction step which fails to yield a feasible control- 
parameter vector. 

Finally the algorithm has two stopping conditions. First the search is 
stopped if both the change in the cost function and the length of the change 
in the control-parameter vector between two successive optimization steps fall 
below their respective input tolerances. Symbolically 



and 


where 11 and ij are the control-parameter vectors resulting from the optimiza- 
tion steps in two consecutive iteration pairs. Second, the procedure is termi- 
nated if the maximum permissible number of iterations specified by the user is 
exceeded. 

The decomposition approach is based on partitioning of the total mission 
to an ordered sequence of self-contained mission segments such as ascent, re- 
entry, and so on. The constraints in each mission segment are then used to 


v+2 v < e 

u - u u 
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define an associated full rank targeting subproblem. Sequential solution of the 
subprob lems ensures that the majority of the mission constraints are satisfied 
on each master level iteration. A master problem, which represents the complete 
problem, is then formulated in terms of the individual subproblem performance 
functionals, constraints, and control variables. Included in the master prob- 
lem are intersegment constraints that cannot be satisfied at the subproblem 
level. (Total AV budgets or propellant limits are .examples of master-level 
constraints that couple subproblems.) The master problem is optimized using a 
two-level procedure with the master-level algorithm optimally coordinating 
solution of the subproblems. The dual-level algorithm implemented employs the 
accelerated gradient-projection algorithm to solve the nonlinear program de- 
fined by the master problem, and the Newton-Raphson algorithm to solve the sets 
of nonlinear equations defined by the subproblems. 

A unique consequence of using a two-level algorithm in conjunction with 
this decomposition approach is that intermediate target values, which are held 
fixed at the subproblem level, can be used as independent variables to be op- 
timized as the master level. 


PROGRAM USE 


This section provides the user with the basic procedures required to use 
the program. At this point, it is assumed that the user has studied the basic 
program capabilities and has concluded that the program can be used to solve 
his problem. The question to be answered next is what does the user have to 
do to execute the program. 


Required Preliminary Analysis 

The user must first define the problem as a sequence of events. This can 
best be done by constructing a worksheet containing the event descriptions in 
columnar form. The user should then assign an integer to each event beginning 
with the first event, say 10, and incrementing by an arbitrary integer, say 
10, for subsequent events. There must be at least two events for each prob- 
lem, but there is no upper limit to the number of events. 

The user must specify the condition at which each event is to occur. This 
should be done by writing the name of the condition and the desired value beside 
each event. The first event begins at the initial conditions specified by 
input and does not require this information. 

The user should next identify the vehicle characteristics and specific pro- 
gram options required at each event-. Any required tabular data such as aero- 
dynamic coefficients, etc, should be included. This identification can be the 
users own terminology. 


53 



The user can now translate the worksheet to program input with a minimum 
of effort by proceeding according to the discussions that follow. A typical 
schedule associated with this translation is given in Figure 21. 


POST Input Phases 

An execution of the program consists of (1) processing the input data for 
all phases, (2) checking for input errors, (3) initializing the equations of 
motion, (4) propagating the trajectory until interrupted by the occurrence of 
the user-specified condition for the next event, (5) reinitializing the equa- 
tions of motion with new user inputs for the event causing the interruption, 

(6) repeating steps 4 and 5 until the user-specified final event is reached, 

(7) terminating the problem, and (8) processing subsequent cases, if any. A 
phase is defined as the time increment from the occurrence of an event to the 
occurrence of the next event. As a result, the names phase and event are 
used interchangeably in the following discussions. 

The first step is to construct a skeleton deck of namelists in the proper 
sequence. Namelist $SEARCH is required at the beginning of each problem. 

This is followed by one namelist $GENDAT for each event defined on the work- 
sheet. If a given phase requires table input, one namelist $TBLMLT must 
follow the $GENDAT namelist for that phase. The $TBLMLT namelist is followed, 
in turn, by one namelist $TAB for each table to be input. This basic NAMELIST 
sequence, which must be followed, is illustrated in Figure 22. This completes 
the skeleton deck setup. The procedure now is to translate the data described 
in the worksheet to program input variables. 

The first step in supply the inputs is to translate the event conditions 
to program inputs. Three namelist $GENDAT inputs are required to accomplish 
this task for each event. First, the input variable EVENT is set equal to the 
event number listed on the worksheet (e.g., EVENT » 1.0). Second the Hollerith 
input variable CRITR is required for all events except the first and is set 
equal to the program symbol for the desired condition (e.g., CRITR = 6HALTIT0, 
for altitude above the oblate planet) . Third, the input variable VALUE is re- 
quired for all events except the first and is set equal to the desired value 
at which the event is to occur (e.g., VALUE = 10 000, for a value of 10 000 ft 
or meters depending on whether English or metric units were selected in namelist 
$SEARCH by setting the variable I0FLAG to the desired value) . On completion 
of these inputs, the user has only to specify the vehicle characteristics, the 
steering options and the desired program options before the job is ready to be 
executed. A summary of key POST input rules is shown in Figure 23. 
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Figure 21.- POST Data Deck Setup Logic Flow Diagram 
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Figure 22.- POST Input Structure 
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Figure 23.- Summary of Key POST Input Rules 
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Namelist Features 


All program inputs are made using the namelist input feature. 

The namelist input format allows card columns 2 through 80 to be used for 
input. Each variable input must be followed by a comma before the next vari- 
able is input. The variables are input as follows: 

NAME1 * VALUE 1 , NAME 2 = VALUE2, ETC. 

Variables may be input one to a card or several to a card, depending on user 
preference. 

Subscripted variables may be input as an array or as individual elements. 
The first element of the array is assumed if no subscript is present on a sub- 
scripted input variable. The following example shows various methods of input 
that are all equivalent: 

EVENT => 1.0, 1.0, 

or 

EVENT (1) = 1.0, EVENT (2) * 1.0, 

or 

EVENT = 1.0, EVENT (2) * 1.0. 

In general, decimal points should not be used with integer type variables. 
In any case, if no decimal, point is input for a variable, the decimal is as- 
sumed to be after the last digit for that variable. As a result, it is best to 
omit decimal points for all variables unless required. 

The namelist used in POST is an extension of the standard FORTRAN namelist 
and has the following added capabilities. 

!) Special List Options .- The initial dollar ($) for each namelist input 
may be preceded in column 1 by any of three print control characters 
as follows: Blank will produce a card image print only for cards on 

which errors were detected; P will produce a print of each card en- 
countered on the input file; and L will print each card encountered 
on the input file, but will also insert the card count after the card 
image . 

2) Embedded Comment Cards . — C or slash in column 1 of any card will cause 
the card to be treated as a comment. Furthermore, when a slash is en- 
countered in any other column, the slash and the remainder of the card 
to the right of the slash will be treated as a comment. 
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3) Special Hollerith Strings .- Hollerith strings may be input as 
OHXSTRINGX where STRING is the Hollerith character string and X is 
any character not contained in STRING. For example, to input the 
Hollerith characters ABCDEF111X in variable W, W = 0H*ABCDEF111X* , 
would be acceptable. The symbol 0 for this option is the number zero. 

4) Repeated Specification .- The following notation may be used to repeat 
a parameter value in any number of successive array locations! M*XXX, 
where XXX can be either REAL, INTEGER, HOLLERITH, LOGICAL, or OCTAL 
input parameter values and where M is the number of successive loca- 
tions in any array where XXX is to be stored. For example, X = 10*1., 
would set X(l) through X(10) - 1. 

5) Improved Error Diagnostics .- The extended namelist prints a diagnostic 
below any card in error with an arrow pointing to the erroneous column 
on the card . 


General Data Input (GENDAT) 

The general data inputs are the constant valued variables and arrays that 
can be input in any phase. These exclude table data and table multipliers that 
are described later. All general data inputs are made via namelist ($ GENDAT ) . 

The general data inputs include the following categories: 

1) Aerodynamic Inputs 

2) Aeroheating Calculation Inputs 

3) Atmosphere Model/Winds Inputs 

A) Event Criteria/Phase Definition Inputs 

5) Gravitational Inputs 

6) Methods of Guidance (Steering) Inputs 

7) Initial Position and Velocity Inputs 

8) Numerical Integration Method Inputs 

9) Program Control (NPC) Inputs 

10) Propulsion/Throttling Inputs 

11) Vehicle/Propellant Weight Inputs 

The selection of program options is accomplished by the user assigning par- 
ticular values to the program control array, NPC(I). A summary of the program 
options controlled by this array is given in Table 11. The elements of this 
array are used internally to select the appropriate calculations for the 
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TABLE 11.- SUMMARY OF THE PROGRAM CONTROL FLAGS 


13 

Units 

Stored 

Value 

Definition 

mm 

Integer 

0 

Conic Calculation Flag 


Integer 

1 

Integration Method Flag 

w£j. iLijjj 

Integer 

A 

Velocity Vector Initialization Flag 


Integer 

2 

Position Vector Initialization Flag 

| 

Integer 

2 

Atmosphere Model Flag 

NPC (6) 

Integer 

0 

Atmospheric Winds Flag 

NPC (7) 

Integer 

0 

Acceleration Limit Option Flag 

NPC (8) 

Integer 

1 

Aerodynamic Coefficient Type Flag 

NPC (9) 

Integer 

0 

Propulsion Type Flag 

NPC (10) 

Integer 

0 

Static Trim Option Flag 

NPC (11) 

Integer 

0 

Functional Inequality Constraints Option 
Flag 

NPC(12) 

Integer 

0 

Crossrange and Downrange Option Flag 

NPC (13) 

Integer 

0 

Propellant Jettison Option Flag 

NPC(14) 

Integer 

0 

Hold-down Option Flag 

NPC(15) 

Integer 

0 

Aeroheating Rate Option Flag 

NPC (16) 

Integer 

0 

Gravity Model Option Flag 

NPC (17) 

Integer 

0 

Weight Jettison Option Flag, Based FMASST 

NPC (18) 

Integer 

0 

Trajectory Termination Flag 

NPC (19) 

Integer 

1 

Flag to Control Printing of Input Con- 
ditions for Each Phase 

NPC (20) 

Integer 

0 

Flag to Specify the Type of Special In- 
tegration Step Size (DT) Prediction to 
be. Used 

NPC (21) 

Integer 

0 

Flag to Indicate the Method by which Flow- 
rate is to be Computed for Rocket Engines 

NPC (22) 

Integer 

0 

Throttling Parameter Input Option Flag 

NPC (23) 

Integer 

0 

Flag that Controls the Velocity Margin 
Calculations 

NPC (24) 

Integer 

0 

General Integration Variable Flag 

NPC(25) 

Integer 

0 

Velocity Loss Calculation Flag 

NPC (26) 

Integer 

0 

Special Aeroheating Calculations Flag 

NPC (27) 

Integer 

0 

Activation Flag for the Option to In- 
tegrate the Flowrate of Selected Engines 

NPC(28) 

Integer 

0 

Tracking Station Option Flag 

NPC(29) 

Integer 

0 

Analytical Vacuum Impact Point Calculation 
Flag 

NPC (30) 

Integer 

0 

Weight Calculation Option Flag 

NPC (31) 

Integer 

0 

A Flag to Activate the Vernal Equinox, Sun- 
Shadow, and Sun Angle Option 

NPC (32) 

Integer 

0 

The Parachute Drag Option Flag 
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options selected. Thus, although the program contains numerous options, any 
user who does not intend to use all program capabilities is not penalized in 
terms of run time. The inputs corresponding to the options selected via the 
NPC array must be specified as input in namelist GENDAT unless the stored 
values are to be used, in which case the variables need not be input. An ex- 
ample of a typical GENDAT setup is given in Figure 24. 

The selection of the guidance and steering options is accomplished by the 
user assigning particular values to the guidance control array, IGUXD(I) . The 
general structure of this array is depicted in Figure 25, and Table 12 con- 
tains a summary of these guidance and steering flags. Admittedly, the hier- 
archical properties of the IGUID array makes -it more difficult to learn than 
the NPC array. The reason for this additional complexity is that guidance 
flags are interrelated. This means that having selected a value for one IGUID 
implies that values must also be assigned to certain other specific IGUID ele- 
ments. Over the years, other approaches have been evaluated, and this decision 
tree approach seems to be the most efficient way to provide flexibility with 
the minimum number of program control flags. 


Tabular Data Input (TAB and TBLMLT) 


The tabular data inputs are the tables and table multipliers that can be 
input in each phase. 

Each table has its own numeric multiplier. All numeric table multipliers 
are input in namelist TBLMLT. The input symbol for a numeric multiplier is 
always formed by replacing the T at the end of the table name by an M which 
denotes multiplier. For example, the multiplier for table PREST is PRESM. All 
numeric multipliers are preset to 1.0 in the program unless overriden by input. 


In addition to the numeric table multipliers, there are some special pur- 
pose multipliers for certain aerodynamic tables that are specified by Hollerith 
names. These are called MNEMONIC multipliers and can be any internally com- 
puted variable that is in the output variable list. These multipliers are used 
to multiply table lookup values of the corresponding aerodynamic coefficient 
tables. All MNEMONIC multipliers are input in namelist TBLMLT. The list of 
tables that have this feature and the corresponding MNEMONIC multipliers are 
as follows: 


Table 

name 

Mnemonic multiplier 
input symbol 

Stored 

value 

Probable value if constant valued or 
monovariant tables are input 

CADPT 

CADPNM 

DFLP 

DFLP 

cadyt 

CADYNM 

DFLY 

DFLY 

CAT 

CANM 

OSE 

ALPHA OR ALPTOT 

CDDPT 

CDDPNM 

DFLP 

DFLP 

CDDYT 

CDDYNM 

DFLY 

DFLY 

CDT 

CDNM 

ONE 

ALPHA 

CLDPT 

CLDPNM 

DFLP 

DFLP 

CLT 

Cl.NM 

<TNF. 

ALPHA 

CMAT 

CMANM 

ONE 

ALPHA 

CMDPT 

CMDPNM 

DFLP 

DFLP 

CNAT 

CNANM 

ONE 

ALPHA 

CNDPT 

DNCPNM 

DFLP 

DFLP 

CWBT 

CWBNM 

ONE 

BETA 

CWDYT 

CWDYNM 

DFLY 

DFLY 

CYBT 

CYBNM 

ONE 

BETA 

CYDYT 

cydynm 

DFLY 

DFLY 
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Figure 24.- Summary of General Data Inputs 


Attitude Calculation 
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Figure 25.- IGUID Array Structure 











TABLE 12.- SUMMARY OF THE GUIDANCE AND STEERING FLAGS 


Input 

Symbol 

Units 

Stored 

Value 

Definition 

IGUID(l) 

Integer 

0 

Type of Guidance (Steering) 
Desired , i.e.. Body Rates, 
Aerodynamic Angles, or Euler 
Angles 

IGUID (2) 

Integer 

0 

Attitude Channel Selector, 

IGUID (3) 

Integer 

0 

Flag to Specify the Steering Option 
when Commanding All Channels Simultaneous- 
ly Using Aerodynamic Angle of Attack, 
Sideslip, and Bank 

IGUID (4) 

Integer 

0 

Euler Engle Steering (Inertial or Relative) 

IGUID (5) 

Integer 

1 

Aerodynamic Angle Rate/Inertial Body Rate 
Combinations Flag 

IGUID (6) 

Integer 

0 

Separate Channel Options for Angle of 
Attack 

IGUID(7) 

Integer 

0 

Separate Channel Options for Side Slip 
Angle 

IGUID (8) 

Integer 

0 

Separate Channel Option for Bank Angle 

IGUID(9) 

Integer 

0 

Separate Channel Option for Yaw Angle 

IGUID (10) 

Integer 

0 

Separate Channel Option for Pitch Angle 

IGUID (11) 

Integer 

0 

Separate Channel Option for Roll Angle 

IGUID (12) 

Integer 

2 

Inertial Body Rate Initialization Flag 

IGUID (13) 

Integer 

1 

Yaw Reference Option Flag 

IGUID(14) 

Integer 

0 

General Open or Closed Loop Guidance Option 

IGUID (15) 

Integer 

0 

General Open Loop Guidance Override Option 
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Some examples of how these MNEMONIC table multipliers are used are as 
follows: 

1) Assume that the normal force coefficient slope is given per degree 
ALPHA. In this case, CNAT would be input as the coefficient slope 
with the CNANM input as 

CNANM - 5H ALPHA , 

the CNAT table lookup value would then be multiplied by the value of 
ALPHA to obtain the value of the normal force coefficient. 

2) Assume that the normal force coefficient is given as a bivariant func- 
tion of ALPHA and MACH. In this case, the MNEMONIC multiplier CNANM 
would be input as 

CNANM = 3H0NE , 

the normal force coefficient would then be the actual table lookup 
value . 

Each table is input in namelist TAB as the input array called TABLE. As 
a result, each table being input requires a separate input of namelist TAB. 
Termination of the table input data for a given phase is indicated by the pres- 
ence of ENDPHS=1 in the last input of namelist TAB for that phase. 

The table inputs for POST are generalized to include: 

1) A preset allowable size for each table. The total size of all tables 
is limited only by the core storage allocated for tables. Both of 
these values can be changed by a simple program modification to satisfy 
user requirements. 

2) Generalized argument specification. The argument to be used for each 
table is specified by input and can be any computer output variable. 

3) Constant valued, monovariant, bivariant, or trivariant table types. 

4) Linear or cubic interpolation capability. 

An example of typical $TBLMLT and $TAB input is given in Figure 26. 


Targeting and Optimization Input (SEARCH) 

The targeting and optimization inputs are made via namelist $SEARCH. 

POST has the capability of performing targeting with or without inequality 
constraints, unconstrained optimization, and constrained 9equality) optimiza- 
tion. The generality of POST enables the user to select the independent and 
dependent variables for the problem from a list of over 400 program variables. 
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For an optimization problem, the optimization variable must be specified 
by input. The variables OPTVAR, OPTPH, and OPT are used for this purpose. 

For targeting (constrained) problems, the dependent variables must be speci- 
fied by input. The variables NDEPV, DEPVR, DEPPH, IDEPVR, DEPVAL , and DEPTL 
are used for this purpose. Both sets of inputs are needed for a constrained 
optimization problem. In any case, the search mode (SRCHM) , the number of 
iterations (MAXITR) , and the independent variables (NINDV , INDVR, INDPH, and 
U) must be specified by input. 

In addition to the above required namelist $SEARCH inputs, there are sev- 
eral others that can be used to further increase the rate of convergence on 
difficult optimization problems. For the most part, these inputs are related 
to problem scaling (MODEW and WOPT) , search direction stepsize control (PCTCC) , 
and convergence tolerances (CONEPS(I)). 

Any type of event (primary, secondary, or roving) can be used in targeting/ 
optimization. However, the user must ensure that the events selected will 
always occur. The association of an event number with the definitions of the 
targeting and optimization variables enables such things as intermediate tar- 
geting and optimization to be performed with the program. This correspondence 
also enables the program to remember the state variables at the beginning of 
the phases where the independent variables are introduced. Thus, when integra- 
ting the perturbed trajectories and the trial steps, only that segment of the 
trajectory affected by the control parameters being changed is integrated, 
thereby reducing the time required to generate the sensitivity matrix. 

The program has the capability to either minimize or maximize a specified 
variable with or without satisfying specified constraints. 

The following inputs are required to define the optimization variable. 

All variables associated with the optimization process are input in namelist 
$ SEARCH. 

1) The Hollerith name of the optimization variable. Input as the variable 
OPTVAR. Any of the output variables can be used as an optimization 
variable, provided it is sensitive to at least one of the independent 
variables . 

2) Type of optimization desired - OPT. If no optimization is desired, 
input OPT as 0. If the variable defined by OPTVAR is to be minimized, 
input OPT equals -1. If OPTVAR is to be maximized, input OPT equals 

+ 1 . 

3) The phase number (EVENT) associated with the optimization variable - 
OPTPH. The variable specified by OPTVAR will be optimized at the piu^ 
side of the event specified by OPTPH. 

4) The weighting for the optimization variable, WOPT. WOPT should be in- 
put as approximately 1 over the anticipated value of OPTVAR. 
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5) The decimal percentage change to be made in the control variables for 
the initial trial step on each iteration - PCTCC. 

The dependent variables can be either equality constraints or inequality 
constraints based on user input. Any calculated output or independent vari- 
able can be used as a dependent variable and as many as 25 can be specified. 

To .constrain a control (independent) variable, declare it to be a dependent 
variable and set IDEPVR non-zero. For example, suppose ALPHA is to be a con- 
trol variable in phase N, but must not exceed 45 degrees. Set INDRV(I) equals 
5HALPHA, INDPH(I) equals N, DEPVR(J) equals 5HALPHA, DEPPH(J) equals N, 
DEPVAL(J) equals 45., and IDEPVR(J) equals 1. 

The inequality constraints may be either functional or single-valued. 

All inputs associated with the dependent variables are input in namelist 
$SEARCH. The required variables are as follows: 

1) The Hollerith name of each dependent variable. Input in the array 
DEPVR(I) , 1-1, NDEPV. 

2) The phase number (EVENT) associated with each dependent variable - 
DEPPH(I) , 1=1, NDEPV. The variable specified by DEPVR(I) will be 
satisfied at the plus side of the event specified by DEPPH(I) . If 
DEPPH(I) is greater than FESN or not input, the last phase is assumed. 

3) The desired value of each dependent variable - DEPVAL(I), 1=1, NDEPV. 

4) The desired accuracy tolerance in which DEPVR(I) is considered to be 
satisfied - DEPTL(I) , 1=1, NDEPV. 

5) The number of dependent variables to be used - NDEPV. The first NDEPV 
variables in the array DEPVR(I) will be used as dependent variables. 

6) The type of constraint for each DEPVR(I) - IDEPVR(I), 1=1, NDEPV. 

The following values for IDEPVR(I) can be specified - 

IDEPVR(I) * 0, if DEPVR(I) is to an equality constraint 

“ if DEPVR(I) is to be an inequality constraint with 
an upper bound 

* “ 1» if DEPVR(I) is to an inequality constraint with a 
lower bound. 

The search/optimization option allows the user to specify as many as 25 

variables for each problem. The initial value of each control variable 
and the phase in which it occurs are also specified by input. If a control 
variable is a guidance (steering) variable, such as pitch rate or angle, the 
appropriate guidance option (IGUID) must be requested in namelist GENDAT for 
the corresponding phase. 
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The value of a given control variable as calculated by the targeting/ 
optimization algorithm will be carried over from one phase to the next until 
overriden by user input or by a new control variable. For example, if the 
linear term in the pitch angle polynomial (Hollerith input symbol PITPC2) is 
a control variable in phase 1.0, and is not a control variable in phase 2.0, 
the calculated value of the pitch rate in phase 1.0 will continue into phase 
2.0 unless the coefficient PITPC(2) is input in namelist GENDAT for phase 2.0. 

A control variable may be constrained by also defining it to be a depend- 
ent variable with an upper or lower bound. 

The control parameters can be selected from any variables in the follow- 
ing categories: 

1) Variables in namelist GENDAT such as initial vehicle position and 
velocity, initial vehicle orientation, vehicle attitude polynomial 
coefficients, etc. For example, suppose the initial velocity (VELI) 
is to be used as a control parameter. The inputs would then be as 
follows: 

INDVR (1) - 4HVELI, 

INDPH(l) - 1, 

U(l) “ AA. , 

2) Constant valued table multipliers in namelist TBLMLT. For example, 
suppose the table multiplier for the thrust table (TVC1M) is to be 
used as a control parameter. The inputs would then be as follows 

INDVR (1) - 5HTVC1M, 

INDPH(l) = XX., 

U(l) - AA. , 

3) User-specified Y arguments from user-specified tables in namelist TAB 
the tables and Y values to be used are specified in the variables 
TABL and TABLY which are input in namelist SEARCH. This feature is 
limited to monovariant tables. For example, suppose the first and 
fourth Y arguments of the monovariant PITT table are to be used as 
control parameters. The inputs would then be as follows — 

TABL(l) = 4HPITT, 4HPITT, 

TABLY (1) = 1,4, 

INDVR (1) = 6HTABL1, 6HTABL2 , 

INDPH(l) - XX., YY., 


U(l) = AA., BB., 
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Each control parameter is uniquely specified by the following variables, 
which are input in namelist SEARCH, and are as follows: 

1) The Hollerith name of each control variable. This is input as the 
array INDVR(I) , 1=1, NINDV. 

2) The initial value of each control variable. Input as the array U(I) , 
1=1, NINDV. If SRCHM is non-zero, U(I) overrides any value input for 
that variable in namelist GENDAT for phase INDPH(I) . 

3) The phase number (EVENT) at which each INDVR is initiated. Input as 
the array INDPH(I), 1=1, NINDV. The variable whose Hollerith name 
appears in INDVR(I) is set equal to U(I) at the beginning of the phase 
specified by INDPH(I) . 

4) The perturbation for each control variable to be used to generate the 
sensitivities. Input as the array PERT(I) , 1=1, NINDV. The sensi- 
tivity DE (J) /DU (I) is determined by finite differencing U(I) by PERT(I) 
and calculating the change in each E(J) . For variables whose nominal 
value is greater than 10.0, the value of PERT(I) should be input 
roughly six orders of magnitude less than the nominal value of the 
variable. The stored values for PERT(I) are 1.0E-4. 

5) The number of control (Independent) variables to be used. Input as 
the variable NINDV. The first NINDV variables in the array INDVR(I) 
will be used as control variables. The number of control variables 
must be greater than or equal to the number of target variables plus 
the optimization variable unless some of the target variables are in- 
equality constraints. An example of a typical setup for namelist 
$SEARCH is presented in Figure 27. 
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$ SEARCH Data Usually Required 
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Figure 27.- Optimization Inputs 



SAMPLE CASES 


This section presents six example problems that exemplify 3 DOF ascent, 
entry and orbital optimization problems; as well as, a 6 DOF simulation of 
the Space Shuttle entry. Each sample caae was selected to illustrate a key 
feature of the program. All examples were taken from actual trajectory prob- 
lems and are not especially constructed sample cases. In fact, the data decks 
shown are the actual working data decks with the removal of redundant data cards 
being the only cosmetic clean-up. Hopefully, by studying these example cases, 
a new user will be able to correlate the details contained in references 2 and 
4. In addition, the examples provide general guidelines that can be followed 
when setting up similar problems. 


Example 1. Space Shuttle Ascent Optimization 

An important ascent trajectory optimization problem during the conceptual 
phases of the Space Shuttle program was that of determining the maximum payload 
capability of various configuration concepts. One such Space Shuttle configura- 
tion is shown in Figure 28. The four key components of this configuration are 
the orbiter, two solid rocket boosters (SRBs) , and the external tank (ET) . 

This multibody nonsymmetrical configuration created special simulation require- 
ments that motivated many of the features contained in POST. For example to 
accurately predict the performance capability of an unsymmetrical configura- 
tion, such as Space Shuttle, it is important to include the thrust vectoring 
losses encountered as the engines gimbal to balance the aerodynamic moments 
caused by the configurations nonsymmetrical shape. This fact led to the devel- 
opment of the static trim option employed in this sample case. 

Trajectory Profile and Problem Formulation.- There are a number of ways to 
formulate the problem of maximizing payload for a given configuration. Each 
approach is based on (1) what is known about the configuration, and (2) what is 
known about the basic trajectory to be flown. In this first example, it is 
assumed that the user knows the dry weight and the propellant weight of each of 
the four major components of the vehicle. Assuming that all the propellant is 
consumed during the flight, which is ensured by terminating the simulation on 
the event criteria 


W 

prop 


= 0 , 


enables the payload weight to be computed from the equation 


(35) 


W 


PLD 


W B0 - “dry’ 


(36) 


Where W is the total burnout weight 

dU 

known weight of the remaining vehicle 
a given configuration, maximizing W 

BO 


(at the final event) and W, is the 

dry 

components. Because W, is constant for 

dry 

is equivalent to maximizing W . Thus, 

x L»D 
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in this example the optimization variable was selected to be W , which is 

BO 

computed as the weight of the vehicle at the instant that the weight of 
propellant is zero. 

As in any trajectory problem, there are a variety of ways in which to simu- 
late this mission, and Figure 28 illustrates the approach taken in this example. 
As indicated, the simulation starts with a 15 second vertical rise, followed by 
a sequence of constant pitch rate steering segments. The static trim option 
is used during all early flight phases, and a three-g acceleration limit is 
enforced after 60 seconds of flight. Event 8.0 specifies burnout of the SRBs, 
rhich are jettisoned seven seconds later at Event 9. Notice also that in 
Event 9.0 new configuration data for propulsion and aerodynamics are input. 

‘These data represent the orbiter plus ET combination that are used for the re- 
mainder of the trajectory simulation. As mentioned earlier, the final event 
criteria is the weight of propellant. Because the last initialization of 
weight of propellant was in Event 9.0, the program variable ^ represents 

the amount of propellant in the orbiter plus ET combination at any time. Thus, 
the final condition 


W 

prop 


0 , 


limits the amount of propellant that can be consumed in all flight phases after 
the occurrence of Event 9.0. 


In this example, the mission requirements are the delivery of the payload 
to the perigee of a 50 x 100 n mi parking orbit. These requirement^ are math- 
ematically equivalent to the three terminal equality constraints 

h - 303 805.0 ft 

V - 25 853.0 fps 



where the subscript, f, denotes final burnout conditions. Extensive computa- 
tional experience indicates that the |h, V^, y^ j constraints are easier to 

satisfy than are their orbital counterparts |h^, h^, 0 ) . The reason for this 

is probably related to the nonlinearities involved, with / h, V , y^\ appear- 
ing more linear in the independent variables. \ / 

Finally, the control parameters selected are eight inertial Eulerian pitch- 

rates, and the initial gross weight of the vehicle at lift-off. Six of these 

pitchrates are used to steer the vehicle during the SRB boost phases, and two 

are used during the exoatmospheric phases. The motivation for using these 

particular control variables is computational experience, which shows that pitch 

angle steering is an efficient technique for optimizing ascent trajectories. 

The particular "breaktimes" in the pitch history were selected after a few 

single pass simulations. The initial gross weight of the vehicle, W , is 

G 
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employed as an independent variable to maximize the payload because in this 

set-up there is a direct one-to-one correspondence between an increment in 

W and an increment in W because all vehicle dry weights and propellant 
G rLD 

weights are held constant during the optimization. 

The previous discussion can be summarized by stating the precise mathe- 
matical formulation of the problem: Determine the control parameters 

U « (W G , 0 lf 9 2 9 ^ 3 * ^q), ( 37 ) 


that maximize: W 

BO 

Subject to: - 303 805 ft ■ 0 

V If - 25 853 fps - 0 
Y If - O' - 0 

The input data deck for this sample case is presented in Figure 29. It is 
important to understand the correspondence between the mathematical formulation 
presented and the specific input data contained in $SEARCH in Figure 29. Table 
13 is a summary of the key targeting and optimization variables as a function of 
the iteration number. 

TABLE 13.- ITERATION SUMMARY FOR SPACE SHUTTLE SAMPLE CASE 


Iteration 

Optimization 

Variable, 

“bo ' lb 

Constraint 
Error, ^ P 2 

Optimization 
Indicator, ^ CTHA 

0 

308 

000 

3.98 

X 

10 7 

— 

1 

305 

151 

4.48 

X 

10 3 

— 

2 

304 

882 

6.37 

X 

10- 1 

89.903 

3 

308 

320 

7.02 

X 

10 1 

— 

4 

308 

369 

4.17 

X 

10- 3 

89.964 

5 

309 

328 

2.21 

X 

10" 1 

89.976 

6 

310 

278 

6.09 

X 

10" 1 

89.999 

7 

310 

300 

5.32 

X 

10 -5 

— 
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P* SEARCH 

C ****************************************************************** 
C PROBLEM - MAXIMIZE WEIGHT 


c 


SUBJECT 

TO - ALTITO 

- 303805 

= 0 


c 



VELI 

- 25853 

= 0 


c 



GAMMA I 

0 

= 0 


c ****************************************************************** 

SRCHM 

X 

4, 

IPRO 

-It 

MAXITR 

* 10, 

OPT 

x 

1 , 

OPTVAR = 

6 HWEIGHT 

, OPTPH 

* 12.0, 

WOPT 

x 

l.OE— 6* 

CONEPS 

89.98, 



NINOV 

X 

V , 

PERT * 

1 .0, 



INOVR 

S 

6 HWGTSG , 

6HPITPC2, 

6HPITPC2, 

6HPITPC2, 

6HPITPC2, 

INDPH 

X 

1 , 

2 , 

3, 

4, 

5, 

U 

*4031000.0 , 

—1 .8 , 

-.5, 

-.2, 

-.3, 

INDVR 16 ) 

x 

6HPITPC2, 

6HPITPC2, 

6HPITPC2, 

6HPITPC2, 


INDPH (6) 

X 

6 , 

7, 

9, 

10 , 


11(6 ) 

c 

-.25, 

-.3, 

-.15, 

-.05, 


NDEPV 

X 

3, 





DEPVR 

X 

6HALTIT0, 

6 HVELI , 

6 HGAMMAI, 



DEPVAL 

=■ 

303805.0, 

25853.0, 

0 .0, 



DEPTL 

X 

100 .0, 

• It 

.001, 




PSGENDAT EVENT 


TITLE 

*5 OH SAMPLE 

PROBLEM 

FOP ASCENT W/ 

DROP TANK ORBITER 

NPC (2 ) 

X 

1, 4, 2, 

NPCC 8 ) * 

2 

t It NPC (161 * 1 , NPCf 21 ) * 1, 

IGUIDll ) 

X 

It 

IGUIDI4) 

X 

1 1 


MAXTIM 

X 

1000 . 0 , 

ALTMAX 

X 

2000000 . 0 , 

FESN * 12, 

DT 

X 

5.0, 

PINC 

X 

20 . 0 , 


TIME 

X 

0 . 0 , 





GDLAT 

X 

28.5, 

LONG 

X 

279.4, 

AZL = 90.0, 

NENG 

X 

It 

WPROPI 

X 

2249000.0, 

ISPV * 439.0, 

GXP 

X 

218.42, 

GYP 

X 

0., 

GZP = 33.33, 

SREF 

X 

4500.00 t 

LREF 

X 

218.833, 



PATBLMLt TVC1M 
S 

P 8 TAB TABLE 
t 

Pi TAB TABLE 
$ 


1 .* 

= 6 HTVCIT ,0,5472000.0, 
* 6HAE1T ,0,232.5, 


PiTAB TABLE = 

- 20 ., 

6HC0T 

t 2 . 

6 HMACH 

, 6 HALPHA , 

12,5 

tltltlt 

1 , 1 , 1 , 1 , 1 , 

0.0, 1.456, .5, 

1.585, 

.7, 

1.598, 

. 8 , 

1.242, 

1 ., 

3.157, 

1.2,2.996 

1.5, 1.816, 2 . 0 , 
- 5 . . 

1.301, 

3., 

.650, 

5., 

.482, 

7., 

.382, 

10., .396 

0 * 0 » • 263 f «5 v 

.338, 

.7, 

. 110 , 

. 8 , 

.302, 

l.t 

.690, 

1.2, .671 

1*5 f • 563* ?. 0 « 
o. • 

.480, 

3* « 

.383, 

5., 

.256 , 

7., 

. 212 , 

10 ., .210 

V • F 

0.0, .180, .5, 

• 1 8 t 

• *7 1 

. 200 , 

• 8 , 

.251, 

l.t 

.495, 

1.2, .502 

Figure 29.- 

Input Data Deck for 

Space Shuttle Ascent Sample 
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1.5, 

K 

.*85, 

2.0, .456, 

3., 

.391, 

5., 

.272, 

7., .231, 

10., .231 

0.0 • 

, 

.263, 

.5, .338, 

.7, 

.110, 

. 8 , 

.302, 

1., .690, 

1.2, .671 

1.5, 

.563, 

2.0, .480, 

3., 

.383, 

5., 

.256, 

7., .212, 

10., .210 

20. 

0.0, 

*1.656, 

.5, 1.585, 

.7, 

1.598 , 

• 8, 

1.242, 

1., 3.157, 

1.2,2.996 

1.5, 

1.816, 

2.0, 1.301, 

3., 

.850, 

5., 

.482, 

7., .382, 

10., .396 

9 

PSTAB 

TABLE 

= 6HCLT 

,2, 

6HMACH 

t 6HALPHA « 

12,4,1,1,1,1 

,1,1, 1,1, 


- 20 ., 


0.0, -1.010, .5,-1.025, 

.7, 

-.99, 

.8, —.815, 

1.,— 1.08 , 

1. 2,-1. 11, 

1.5, -.895, 2.0, -.788, 

3 . ,■ 

— .635 , 

5., -.480, 

7., —.43 , 

10., -.43, 

0 • v 

0«0* *015 v • 5f *04 9 

.7, 

.01 , 

.8, -.045, 

1. , .08 , 

1.2, .038, 

1.59 —.02 9 2 .O 9 -.1089 

c 

3 ., 

-.145, 

5., -.15, 

7., -.15 , 

10. ,-.15 , 

5.9 

O.O 9 .5459 .59 *759 

.7, 

.53, 

.8, .365, 

1., .69 , 

1.2, .638, 

1.59 .43 9 2* 9 .2429 

3 ., 

.11, 

5.* .025 , 

7., . 00 , 

10., .00 , 

20 • 9 

O.O 9 2.135 9 « 5 t 2.24 9 

.7, 

2.09, 

.8,1.595 , 

1., 2.52 , 

1.2,2.438, 

1.5* 1.78 9 2. 9 1.292* 

3 ., 

.875, 

5., .55 , 

7., .45 , 

10., .45 , 

PSTAB TABLE = 6HCMAT 

,1 »6HMACH 

,12,1,1,1, 



0.0* .019* .7* .0218* 

.9, 

.0302* 

1., .023, 

1.2 ,-.011, 

1.5, -.032 

1.8 *— .0395* 2. *-.0419* 

3.,“ 

.0396* 

5. ,-.0187, 

7. ,-.0082 

, 10 . , 0 .0 , 

PSTAB TABLE = 6HXREFT 

, 1 »6HMACH 

,12,1,1,1, 



0.0*137.86* .7*140.05* 

.9, 

136.77, 

1., 147.71 

, 1.2,145. 

52, 

1.5*144.43* 1.8*141.589 

2. ,138.3, 

3., 131.74 

, 5. ,118. 83, 

7. *109.42* 10. *91.91* 

ENOPHS * 1* 

* 






* 

PSGENDAT EVENT = 

2, 

CRITR 

= 6HTIHE 

• VALUE 

* 15.0, 

IGUID(4) = 0* 

ENOPHS * l t 






PSGENDAT EVENT 

3, 

CRITR 

= 6HTIME 

• VALUE 

= 25.0, 

ENOPHS * 1* 






9 

PSGENDAT EVENT * 

4, 

CRITR 

= 6HTIME 

• VALUE 

M 

* 

O 

• 

O 

* 

ENOPHS * 1* 






9 

PSGENDAT EVENT 

5, 

CRITR 

* 6HTIME 

, VALUE 

O 

• 

O 

O 

II 


npc< 7) * l* 

ASMAX * 3.0, 

ENOPHS * 1 , 

S 


PSGENDAT 

ENOPHS 

EVENT 

Z 

1 , 

6, 

CPITR 

* 6HTIME 

• VALUE 

= 120.0, 

9 

PSGENOAT 

EVENT 

s 

7, 

CRITR 

* 6HTIME 

• VALUE 

* 150.0, 


Figure 29.- Continued 
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OT 

* 10 . 0 , 

ENDPHS 

* 1 , 

S 

P 8 GENDAT 

EVENT * 

NPC (7) 

* 0 , 

NPC<9) 

* 0 , 

WEI CON 

* 0 . 0 , 

ENDPHS 

* 1 , 

* 

PSGENDAT 

EVENT * 

DT 

* 20 . 0 , 

PINC 

* 50.0, 

NPC (7| 

* 1 , 

NPC ( 9 ) 

* 1 , 

WJETT 

= 665000.0, 

WPROPI 

- 809000.0, 

1SPV 

* 459.0, 

GXP 

* 142.0, 

GYP 

* 0 . 0 , 

GZP 

« 25.0, 

SREF 

* 4840.0, 

LREF 

= 135.0, 

$ 

P 6 TBLMLT 

$ 


8, CRITR 


6HWPR0P , VALUE 


0 . 0 , 


9, CRITR 


6 HTDURP , VALUE 


* 7 ., 


PSTAB TABLE 

t 

PSTAB TABLE 

S 

PATAB TABLE 
- 20 ., 


= 6HTVC1T ,0,1431000.0, 

= 6HAE1T ,0,154.34, 

* 6HCDT ,2 *6HMACH ,6HALPHA *12,7,1,1,1,1,1,1,1,1, 


0, .024, 
2,. 116, 
4. , 

• 2 9 #024 f 
2 #49f # 1 9 

• 69*0269 * 89*0289 
Bf#092» 3*9t.082t 

.9, .035, 
40, .03, 

1.3, .093, 

1.5, .122 

0, .024, 
2,. 116, 
0 .* 

• 2 f • 02 A t 

2#48 9 #1 y 

• 69*0269 * 89*0289 
3 9 *092 v 3*9 f *082 t 

.9, .035, 
40, .03, 

1.3, .093, 

1.5, .122 

0, .026, 

• 2 9 #026 t 

• 69*0269 * 89*0249 

. 9 , .036 , 

1*39 *092 9 

1.5, .118, 

2,. 106, 

2*409 *091 

9 39*0829 3*9 V *074 

• 40,. 022 

9 


5. 

0,.042, .2, .042, .6, .04, .8, .04?, .9, .076, 1.3, .124, 1.5* .142, 

2. . 124, 2.48, .098, 3, .088, 3.9, .079, 40, .033, 

10 ., 

O*. 076, .2, .076, .6, .08, .8,.l, .9,. 13, 1.3,. 194, 1.5, .192, 

2.. 165, 2.48* .127, 3, .114, 3.9, .095, 40, .057, 

20 ., 

O ,. 36 , .2, .36, .6, .362, .8, .44, .9, .41, 1.3, .39, 1.5, .36 , 2, .32, 

2.48. . 242, 3 , • 2 24 , 3.9, .216 , 40, .238, 

30., 

0«. 36 , .2,. 36, .6, .36, .8, .44, .9, .41, 1.3, .39, 1.5, .36, 2* .32, 
7*48, ,44, 3, .416, 3.9, .4, 40, .3, 


Figure 29.- Continued 
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s 

PSTAB TABLE * 6MCLT ,2,6HMACH ,6HALPHA *12 #7*1 *1*1 * 1» 1*1, 1* 1 * 
- 20 .* 

0*— .0T* .2*— .OB* .6*— .12* .8*— .12* .9*— .12* 1.3,— .12* 1.5*“. 12* 
2*-. 13* 2.48*-.14, 3,-. 12* 3.9.-.1* 40,-. 14* 

4. * 

0*“.07« .2*— .OB* .6*-.12* .8*“. 12* .9*— .12* 1.3*“. 12* 1.5* .12* 
2*“. 13* 2.48*-. 14, 3 *— . 12 * 3.9, -.1, 40,-. 14, 

0. * 

0*.08* .2, .08* .6, .08, .8*.06* .9, .06* 1.3*. 07, 1.5**04* 2*0.0* 
2.48, -.02* 3, -.03, 3.9, -.04, 40, .03, 

5. , 

0 , . 29 * . 2 , . 29 * .6 * .29* .8*. 28, .9*. 28, 1.3, .3, 1.5*. 24* 2**17* 

2.48.. 12, 3*. 09, 3.9, .08* 40, .21, 

10 .* 

0 * . 5 * • 2* .6 * .6 *.49* .8**48* .9**52, 1.3*. 52* 1.5*. 41* 2**33* 
2.48**25* 3**2* 3.9, .15, 40, .4, 

20., 

0,.94, .?*.94* .6, .92,- .8, .9, .9, .94, 1.3, .89, 1.5, .75 , 2**63, 

2.48. . 51* 3, .43, 3. 9, .39, 40, .76, 

30. • 

O , . 94 * .2**94* .6**92, .8**9* .9**94* 1.3*. 89* 1.5*. 75* 2*. 68* 
2.48**67* 3*. 65* 3. 9*. 62* 40**76* 

$ 

PSTAB TABLE * 6HCMAT *0,0.0* 

$ 

PSTAB TABLE * 6HXCGT *1 ,6HWEIC0N, 5,1*1 *1* 

0,87.64* 202250,93.93, 404500,99.68, 606750,104.04, 809000*104.37, 

$ 

PSTAB TABLE * 6HYCGT ,0,0.0, 

$ 

PSTAB TABLE * 6HZCGT ,1,6HWEIC0N,5, 1,1,1, 

0,31.33, 202250,31.5, 404500,31.75, 606750,32.42, 809000,33.83, 
ENDPHS = 1* 

$ 


PSGENDAT 

ENDPHS 

EVENT 

X 

X 

It 

10. CP IT* 

* 6HTDURP * 

VALUE 

= 100.0 

% 

PSGENDAT 
NPC ( 1 1 
ENDPHS 

EVENT 

X 

X 

X 

2, 

1, 

11, CRITR 

* 6HTDURP * 

VALUE 

* 150.0 

PSGENDAT 
ENDPHS 
ENDPRB 
END JOB 

EVENT 

X 

K 

X 

X 

It 

1, 

It 

12* CRITR 

* 6HWPR0P * 

VALUE 

* 0. » 


t 


Figure 29.- Concluded 
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In this table, the function P 2 is a measure of the weighted error in the 
constraints, and is calculated from the equation 


'(h) B0 - 303 805' 

2 

+ 

’ 

( V l) 

BO - 25 853 ' 

2 

+ 

( y i)bo1 

100 

L 0.1 J 

O.OOlJ 


(38) 


Clearly, in the general case P 2 is the square of the norm of the weighted 
error vector e_, and is computed as 


The condition P 2 < 1 is sufficient to ensure that all constraints are satis- 
fied to within their specified tolerances. 


The other important optimization output variable is the angle, CTHA. 
Mathematically, CTHA is defined as the angle between the unconstrained per- 
formance gradient vector and the projection of this vector to the plane that is 
a tangent constraint manifold. CTHA converges to 90* at the constrained op- 
timum because the performance gradient becomes orthogonal to the constraint 
manifold at that point. This angle is computed from the equation 


CTHA 


cos 


-1 ( * • P M 

\l|lllll p £ll/ 


(39) 


where £ is the unconstrained performance gradient and P is the projection matrix 
to the constraint tangent plane. 


Typical output for this type of problem would be (1) the trajectory, block 
printout for the nominal trajectory, (2) the sequence of iteration summaries, 

(3) the trajectory block printout of the final optimal trajectory, and (4) plots 
of all user-requested variables. These types of output are illustrated in 
Figure 30. The plot capability is not contained in POST but rather is gener- 
ated from a profile type written by POST. The frequency and amount of trajec- 
tory printout can also be controlled by the user. In this example, the only 
printout is at the minus and plus sides of each event. These data give a brief 
but useful summary of the trajectory time histories. 


Example 2. Single-Stage-to-Orbit (SSTO) Entry 

This SSTO entry case, illustrated in Figure 31, features several advanced 
program options, that are particularly useful in entry trajectory analysis. For 
example, the trajectory control scheme used is significantly more complex than 
the simple open loop technique demonstrated in Example 1. In addition, several 
of the more advanced program simulation options are employed. Options include 
roving events, multiple time channels, feedback steering, monitor variables, 
print block changes, special variables, and multiple runs. 
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TRAJECTORY BLOCK PRINTOUT 


ft r 


b< 


H 

M 


SAMPLE PROBLEM FOR ASCENT W/ DROP TANK OAtlTtfc 
INWJT UNITS » ENGLISH, OUTPUT UNITS «■ ENGLISH 


INITIAL CONOIT1CNS 

TIME * 0. TIMEO 

GCLAT *• 2.BSQ000ouE*Gl GOLAT 
OCR AD * 2 .0925791 01 *07 ALTITn 

VELA - 0. GAMNAk 


2.6 5uuOoOOE *01 LONG 
O. 

0. AZVE LA 


2 .799000001*02 LUNG I * 2 . T990uo0oE *02 


0. 


PA OOP AN TERMINATION PARAMETERS 

FESN * 12.000 MAX T 1 

THE LAUNCH PAO INERTIAL 4L) FRAME IS OEF INEO bV 

LA1L * 2. £5 0006001*0 1 LC*NL ‘ 2.799G0D0OE *02 Ail 


attracting PLANET MCOtL 

RE ■ 2.092 5 7A1UE *G7 RP 

J2 ■ 0. J3 


l.UOOuOUOOi*Oi ALIM1N * -5 .OOOOOOOOE *u 3 AlTMAX * 2 . uouOuuuoE *06 

9.000000001*01 

2 .092579101*07 OMEGA * 7.292 1 lOOOE-Ot MU * I .90765 39 g( * I o 


G. 


J9 


0. 


SAMPLE PROBLEM FOR ASCENT w/ DROP TANK. OREITEP 
BEGIN PHASE 1.000 

1962 u.s. standard atmgsphepe mooel 

AEftOOYNAMIC COEFFICIENTS SPECIFIED oY DRAG AND LIFT CUE Ff IC IE NTS 
SkE F » * . SOOoOGoOE *03 LREF * 2 • 16t33000 E *02 LREFV 

PROPULSION CALCULATED FOR 1 ROCKET ENGINES 

VEHICLE HEIGHT PARAMETERS 

WGTSG » A .03328 1 1 2E *06 WPLO » 0. HPROP 

GO * 3.21 7AOOOOE *01 

NUMBER Of INTEGRALS FOR THIS PHASE * b 
INTEGRATION SCHEME » FOURTH ORDER MUNGl-KUTTA 

OT * S.OOOOOOOOE *0O P1NC * 2.000000001*01 


2 • 29900000E *06 WJCTT 


USE INERTIAL ROLL 
ROLPC * 
PITPC « 
YANPC * 

Y AW I * 

PIT I • 
ROLI 


, YAW, 

0. 

0. 

0. 

0. 

0. 

0. 


ANO PITCH COMMANDS 


VAWR 

PITR 

ROLk 


QYAW 
DPITCH i 
ORCLL 


THE NEXT EVENT TO OCCUR WILL bt ONE OF THE FOLLOWING 
. 2.000 TVPt • PR INKY 

TOL • 1. OOOOOOOOE -06 MOL * 1 


TAfcLIS ANO MULTIPLIERS FOR THIS PHASE 
CCT » l .OOOOOOUOE *00 LLT 
XREFT t 1 .OOOOOoOOE *00 


DESN 

OESN 

DESN 


VALUE > 1.500000004*01 


.000000001*00 CHAT , 1 . OOOOOOOOt *O0 TVC1T , I .OuOuOOouE *wO AE 1 T 


1 .uuOuouIhjE ♦ uu 


PROGRAM CONTROL FLAGS 


NPC 

NPC 

NPC 

NPC 

NPC 


1 ) 
I 0) 
I 15 ) 
122 ) 
(29 1 


NPC 
2 NPC 
0 NPC 

0 NPC 
G NHL 


I 21 
I 91 
4 16) 
123 ) 
130) 


1 NPC 
1 NPC 
1 NPC 
U NPC 
0 NPC 


4 3) * 
1101 * 
4 17) « 
4 29) * 
431) * 


A NPC 
0 NPC 
0 NPC 
0 NPC 
0 NPC 


( A) 

411 » 
4 10 ) 


GU10ANCE CONTROL FLAGS 

IGUID 4 1 1 •» 1 ICC ID ( 2 I 
IGUID ( 6) * 0 IGUID I 9) 
IGUID 415) * 0 IGUID llo) 

IGULt (22) * G IGUID 1231 


IGUID 4 3) 
IGUID (10) 
IGUID 417) 
IGUID (29) 


U IGUID 4 <•> * 
u IGU 10 4111 « 
0 IGUID (lb) « 
0 IGUID 1 25 ) « 


2 NPC 
0 NPC 
0 NPC 
0 NPC 
0 NPC 


1 IGUID I 31 * 
0 IGUID Cl 2 I * 
0 IGUID (19) « 


( 3) 


x 

NPC 

4 6)* 

0 

NPC 

(12! 


0 

NPC 

(13) * 

0 

NPC 

119) 


1 

NPC 

120) * 

0 

NPC 

(26) 


D 

NPC 

(27) * 

0 

NPC 

(33) 


0 

w>c 

4 .»9 ) * 

0 

NPC 


4 

4 L * 


1 IGUID 4 b) * 

2 IGUID (13) * 
b IGUID 42u) - 


7) 
*►> 
( 21 ) 
(2b) 
435) 


IGUID 1 7) * 
IGUID OA) » 
IGUID (21) « 


SAMPLE PROBLEM FOR ASCENT W/ DROP TANK ORhlTEk 


PROBLEM NC. 


**• PHASE 


TIME 

G. 

TIMES 

ALTITC- 

■9. To 8 37 1 38 f — 07 

glkad 

VELI 

1.39101 1601*03 

GAMMA] 

V£lR 

0. 

gammar 

VI CA 

0. 

GAMMAA 

GAMAC 

c. 

A2VA0 

TMkUST 

A ■ VI 997969E *06 

HE 1GHT 

eta 

l .OOOOOOOOE *L»0 

ETAL 

FT XB 

9.9799T969E *06 

F AX b “ 

FTYB 

0. 

F AVP 

FT 20 

0. 

FA2B - 

CA 

1 . BOGOOOOOE -01 

CD 

CN 

1 .500000001-02 

CL 

CY 

0. 

HE AIRT 

OYNP 

c. 

MACH 

TIME 

i . souoooout *oi 

TIMES 

altitd 

9.33310129E *02 

GCRAO 

VE U 

) .39603 lOOt *v3 

GAMMA I 

VELR 

1.29291 7591*02 

GAMMAR 

VI LA 

1.29291 759E *u2 

GAMMAA 

GAMAO 

-5.723115991-03 

AZVAO 

THRUST 

9.99639b 301*06 

WE 1GHT 

eta 

I .OOOOOOOOt *oo 

ETAL 

FT XF 

9 . 996398 36 £ *06 

F AX f • 

FTYB 

0. 

F AYB 

F T 20 

G. 

F A2b - 

CA 

1.630259A7E-O1 

CO 

CN 

1.262960371-02 

CL 

CY 

0. 

heatrt 

DYNP 

1 .93297183E *01 

MACH 


u. TDURP < 

2 .09237A luE *07 GOLAT 
0. A2VEC1 

0. A2VLLF 

9 . DODOOOoOE *01 A2VELA 
j. DWNRNC 

06 WOOT 

00 1PNULI 
AX 6 
AYb 
AZb 

01 DRAG 
-02 L1F1 

TLHE Al 
RE YNO 


a. 03 32b! 12E • 
l.OOOOOuObt* 


1 . 600000 UOE ~ 
1 .500000001- 


1 .SOOOOOuOE *ol 
2 .09 2c>6?A Jt *07 
3.3O039lulE *00 
6.9624,75971*01 
0.56267597E *01 
2.0061792 II *Ot’ 
3.8A631 u 7AE *t>6 
1.000000001*00 
-1.592u235?E*o9 
0. 

■I . 09b 5 7o06E *03 
1.020919621 -01 
1.3 lu5Au?tol -02 
0. 

1 . 1 6 1 79b 3 1 E -O ) 


TDURP 

GDI A I 

AiVEL I 

AiVELF 

A 2 VELA 

DWNRNC 

WOOT 

1PNULL 

AXB 

AYB 

Aib 

DRAG 

LIFT 

TLHE AT 

RE YNO 


M PHASE l. 

. 300OD00OE *01 
. S A9990 Sa E *01 
.002990071 *01 
.1336906AE *02 
. I3369u6AE *02 

.2A6A6V251 *0A 

* 1 66077734 *01 


,000 

01 NS 2.37690*971-03 

Pk£ s 

2. 1162 1660b *03 

ATEN 

5 . 1 bo TwoOOL *02 

GCLAT 2 .BSoOOOUilE *01 

LONG 

2. 799ouOuoE *02 

L0NG1 

2 . 79>>^t<uu( *u2 

XI 3*00359 bOOE *06 

VXI 

1.323009601*03 

AX 1 

1 .Ob 7b 95 ft9t *00 

YI -1.119296271*07 

VY1 

2.1 902 2029E *02 

ATI 

-6 . 5 71 «»© 3301 *00 

21 9 .9699006 31* 0* 

V21 

13. 

A2 1 

3.6l6>?o23E *00 

CRRNG 0. 

DPR NG 1 

w. 

DP PNG 2 

o . 

WE ICON 1 .990116121 -OB 

WPRQP 

2 .2 AVuoOOOE *06 

ASMG 

1 .23*72 lobi »O0 

IV NULL 0. 

1NCPCH 

0. 

INC YAW 

0. 

ALPHA 0. 

alpdct 

0. 

AL PICT 

0. 

FETA 0. 

bl TOUT 

o. 

OALPHA 

V. 

bNKANG 9.000000001*01 

BNKOCT 

0. 

UAL TLT 

O . 

ROLI 0. 

YAW* 

0. 

ROLbO 

U . 

YAW1 U. 

PITR 

E • 999996 92 f *Ol 

PlTfcD 

u. 

PITI 0. 

F CL R 

u. 

YAWbD 

u. 

AS X I 5.7u2uoB52E*uO 

ASTI 

- j. 99930*. I 31 *ul 

AS 2 i 

1 .bvSSS 7out *01 

.000 ••• 

DENS 2 .312670091-03 

PRES 

2.095a 1391E *03 

A 1 E M 

5.1539ls l 5t*u2 

GCLAT 2.699998561 *Ol 

LCNG 

? . 79 39991*31 *02 

LONUl 

2 . 79^o^669k *u2 

XI 3.023529331*06 

VXI 

1.391 10399) *03 

AX 1 

1 .32o*.Ov9Vt*00 

YJ -1 . fc l 90976*1 *07 

VY1 

1.0799b292t *02 

AVI 

-6 ,25 6U«.Yu9E *S»0 

21 9 ■ 9b 5 39 1 36 1 *06 

Vil 

6. 1076 53914 *01 

A2 1 

9.5911 76351*00 


. lb9*»27SlE-U3 
•591732*04 *0A 
. 1399372llE*o3 


ROLI 
V AH 1 
PI T l 

7S96Ao33L *0* AS X I 


CRRNG 0. 0PRNG1 

HE ICON 1.869703171*03 T^RClP 
1 YNUL L 0. INCPUH 

ALPHA -1.909630811-01 ALPOOT 
BETA 3.101953001-01 bl TOOT 
bNKANG 1.233397161*02 INK DOT 


u. 

2 .06202 9611 *06 ASMG 1.2990 54 53* *00 

u. INC YAM 0. 

ALPTuT 3.99iM*j7#bE -Ol 
UAL Pha-2 .B79^1»».2E *0O 
OALTOT 6 .65<3o56 16E *00 


VAHk 
PI Tk 
ROl R 

.970657921*00 ASVI 


.700299011*02 
(.999992971 *ol 
.bDOOOOOOE *02 
).oI22o»22E*ui 


ROLbO 
PI TbD 
YAHKO 
AS 2 1 


] .9M7R0O99E *01 


END OF PHASE 


Figure 30.- Typical POST Output 
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&A*Plt PRC6LE N F CR ASCENT W/ DRCP 1 ANA CREITER 



••• INPUT SUNNARY f GR PAOLLEM NC . 

SI ARCH/OPT IHUAT ION INPUTS 

NINOV « 9 N 

I NOV* . NOT St . 


WGTSt , 

P1TPC2 , 

l.OOO* 
7.000. 

4.u35?6tl2E«G6 
-6.b5?O045 It -02 
2 .4 0077*001 -07 
3.33333333E *00 
2.40G77*oJt-o7 
3.3>3333>3E -04 
ACT 1 TO . 

000 . 000 . 
3.O30O3OOO£*G3 
1.00000 Out t *02 
0 . 


NOE P V • 3 

P1TPC2 , 

P1TPC2 , 

2 . 000 * 

9.000* 

.-4.o2959nut-vl 

*-1.3o635729|-ol t 

» b.»omM7t«00 l 
I i.33553336fc-05, 
* 6.666666671-04, 
Vt II . 

voo.ooo, 

> 2. 363300001 *04, 

. 1 .OOooOOOOE— 01 i 


SRCHM 

PITPC2 , 

PITPC2 , 

3.000. 

10 . 000 . 

>-*.330SJ62uE -01 * 
-1.16711 7751 -01 , 
.2 .UOOOOUUOt *00* 

2 .oooooooot *ui* 

' /.OuOOOOOOE -O-, 
2.000000001 -03. 
CANNAI , 

900.000. 

u. , 

1.00000000E-03. 

0 . 

0 . 


NOOtw « 
PI7PC2 


-1 .677669b Jt -ol . -0. 7 72* j.’S It - 01 , -..o76 72429i -ol , 
3.COOOOOOG1 *t/G, 3 • 3> 3 3333 21 *Lo* * .oOOGOoOol *0O * 
3.000G0O00t-O4, j .3 33333331 -O*. * .0000000 ot -o* , 


1 . 000000001*00 
I .GOOooOUOC -06 
1 .000000001 -03 
1.000000001-01 
1.000000001-01 


NAX1TA - 
P2N1N « 
CONSX? ■ 
PCTCC *■ 
CCNLP3 * 
GANAX ■ 


0 IRPQ » -1 

1 . JOOOOOOOt -01 STNNP2 * 1 . OGGGOOGOt -01 

1 .tOOOOOOGE *02 flTtRl ■ 1 .OGuGGGOOt -04 

1 .00000 Oool-u* CONE PI » 0.94*troouGf *ol 
1 ■ OOOOOCOOE —0 1 CCNEP> «. 1 .00*0 Gw 001 -O 1 



••• ITERATION NUNbih 1 


CP/I TR 4.03460 750 E *00 
PERT 1 . QovutruOOl *00 

1.000000001 -04 

part i als ha aliito with 

SNAT -2.303660201*00 
1.49390 101 E *06 
PARTIAL* Of Vt L I WITH 

SNAT -3.6S04J 7731 -02 
-3.91391 00 1 E ♦ O 3 

partial s of gannai with 

SNAT -3. 630 69 70 SE -OS 
2.72 T 3 5672 E *01 
till) -*.031000001*00 
1.59072 1 16E -1 1 
G1HAG 4.031000001*00 
02 111 1.491342251*10 

1 *07427t731 *00 
02 HAG 1 .493 1 1 04 7 E * 10 
P0K1 1-1.IT325747E-U3 
-9. 69009096c -04 
PGlNAO 6.07707 1001 -03 
WVEC 2 . 400774001 -o 7 

3.333333331*00 
NAC 3 

1AI 1 

OPIOS 1.953995301-02 
DP20S -5.461771261*00 
STPNaX 1 .000000001 *10 
UN AG 3. 1 1 006422 £ *00 
OUNAG I. 450179761-C1 
Gamas T o. 

P1TRY -3.060000001-01 
P2IRY 3.90 |9066OF*O7 
VPREO O. 

INOVR W&TSt 

PITPC7 
INOfH 1.000 

T .000 

UN) 4.0701 5u07t *06 
-3.195777141-01 
WtlU -2.746602731 *01 
BE PV* ALT I TO 

OE PPH 900.000 

till -2.746607731*03 
PI -3-051 500211-01 

OP WAR WEIGHT 

OP 7 PM 12.000 

CPWAL 3 • 05 1 50071 1*05 
P7 4.479900201*03 


I.OOOOGOOOt -04 
1 .OOGUUOOGE —04 
RESPICT Tt Util 

6.541431*51*05 
2.671955431*06 
RESPECT TC U(ll 

-1.303032671*03 
-7.40352320E *03 
RESPECT TO Ulll 

0.690610 701 *uu 
o.2 3*3u60ol *01 
-3.197442311-1] 
0 . 

4.34262 7t 31 *06 
1 .054719001 *00 

-4. 27020 3* 7E— 03 
-0.329069731-0* 

5.555353361-01 

6.666*66671*00 


1.43012 9761-0 
-3.05150*211-0 
4.479961201 *o; 
o . 

P1TPL2 

pitpc; 

2 .000 

9.000 

-1 .0437673*1 .01 
- 1 . 66 30 40 S 6 t - t 1 
—5 .6394V3 b*l *01 
Vlll 

voo.ooo 

-3.639*93641 *oi 


1 .OOOuOOOOl -04 
1 . OOOUOCMJOE -04 

9.620200311*05 
I. 190504 U £ *06 

-2.4346 34701 *03 
-4.929770701 *03 

1 . .06642301 *01 
5. 103679201 *01 
1. 7763 36641 -u 
0.611 704201 -13 

1.0 337064*91 *00 
* • 20943*3*1 *07 

-3.406460161 -04 
1 . 0 1 904 7 U 0 1 — u 3 

2 . 000000001*00 
2 * 000000001 *01 


1.437965731-01 

0 . 

0 . 

4. *79464311 *03 
PITPC2 
P 1 TfC 2 

3.000 

10.000 

-3.03604901 l -01 
-5.2274T4UI -02 
2. .>>4944421 *01 

GANNA 1 

voo.ooo 

2 . 334944*2 1 -02 


I .000000001-04 

1.236030961*06 

-3.320709501*01 

1 . o * 21 l 2261*01 

7.105427361-12 

9.917442 771*07 

*.033329741-04 

5.000000001*00 

0. 

0. 

0 . 

o . 

PI TPC2 

4.000 

-2.013167941-01 


1 .300000001 -04 

3.613150031*06 

- I .007279331 *o4 

*••* 12 396 Cf *ol 
1.39072U6E-U 

*.4*0607951*00 

3 .0060 94000-03 

3.333333331*00 

0 . 

G . 

6 ■ 

o . 

P 1 7PL2 

5 . oO«j 

->.l 1.522001 -01 


►•LLLIN NL. 1 

1 .OOoOuGGOE— 04 

1.3 >0*3 1301 *06 

-9.19333*761*03 

2 . 37ol6> *li*ol 

1.7763566*1-11 

I .3912796 71 *06 

9.20137 301 L —04 

4 . ooooowoOl *O0 

o. 

0 * 

u . 

o . 

p, lPv.2 

O.U0O 

-..57>3l6?7t -u| 


Figure 30.- Concluded 
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-76.6m (251.3 ft) 
1.0 - Start Entry <? 1219.2 k a (400 000 ft) 



^20.0 Start Pullout @ 990.6 km (325 000 ft) + 0 
'v^30.0 End Rollout after 10 s 

'35.0 Interrupt (3 acc • 0.1 fps 2 

^40.0 Terminal Pullout at dh/dv ■ 5.85 s 

50.0 Terminate Heatrate Control @ 2.2 g 

^ \ 

\ \ 60 . 0 Terminate acc Control 14/s after Initiation 

70.0 Start Final Rollout 
Roving Event on ^C 75 '° B “ k 

Altitude 332.2 km ^_---76.0 Drive Bank Angle to Zero 

(109 000 ft) to Start — 80.0 Roving Event on Mach No. to Change 

New Aeroheating 65.0 Roving*''^ Angle of Attack to 20 deg 
Calculation Event on^ - ^.^ 96,0 Rovln S Event on Altitude 

Flight Azimuth \ to Ter,lru,c * Aeroheating 
Calculation* 

-^lOO.O End Trajectory 

Earth 152 km (50 000 ft)^' <a 152 km (50 000 ft) 


Figure 31.- Trajectory Profile and Sequence of Events for SSTO Entry 
Sample Case 
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The simulation is initiated and entry interface ( 400 000 ft), and termin- 

ated when the vehicle has descended to 50 000 ft. The trajectory control scheme 
attempts to minimize total heat by capitalizing on the ability of the TPS mate- 
rial to withstand a high heat rate for a short period of time. The basic op- 
timization strategy is to attain the prescribed maximum heat rate as early as 
possible on entry into the atmosphere. This limit is then followed until the 
acceleration constraint is encountered. The bank-angle linear-feedback steer- 
ing algorithm is then switched from controlling heat rate to controlling the 
acceleration profile. The acceleration limit is then followed until it is 
necessary to deviate from the limit to achieve the specified crossrange. This 
trajectory scheme is presented in Figure 31. Parametric results show that this 
simple approach provides the minimum total heat for this configuration by mini- 
mizing the entry flight time, subject to the heating rate, the acceleration, 
and the crossrange constraints. 

One of the more difficult problems encountered in this example is that of 
determining the proper bank angle schedule to enable the vehicle to "pull out" 
of its initial entry plunge at the maximum heat rate that can be tolerated by 
the TPS. This "pull-out" maneuver is further constrained because it is im- 
portant to minimize the heat rate "overshoot" when the linear feedback steering 
option is exercised. This is accomplished by using the POST iteration feature 
to solve the pull-out problem, which can be defined as: 

determine the desired bank angle, , that satisfies the heat rate equation 

* (*d) * 4 max’ <«» 

when 

Q = 0. 

In the trajectory set-up, presented in Figure 31, the vehicle is initially 
banked ninety degrees (<j> = 90“) and 4 ^ is the desired bank angle when the 

vehicle reaches the 325 000 ft altitude. As such, 4 > d really represents the 

amount of "un-bank" required to achieve the proper pull-out conditions. 

After the proper pull-out conditions are satisfied, the vehicle descends 
using closed-loop steering to maintain heat rate and limit maximum deceleration. 
This leads to the second problem that must be solved to satisfy the cross- 
range requirements. If the vehicle remains on the maximum deceleration boundary 
too long, it will be unable to attain the required crossrange. If it departs 
too soon, it will be able to satisfy the crossrange, but it will increase the 
flight time, and, hence, the total heat. As a result, the following open loop 
optimization problem is formulated at the "bottom" of the trajectory: deter- 

mine the acceleration boundary departure time, T Q , and the bank angle schedule 

characterized by the parameters {c^ <j >2 413 4 ) 4 } 

that maximize: T^ (i.e., stay on g-boundary as long as possible) 
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subject to: C R 1100 n mi = 0 

at the trajectory termination condition h = 50 000 ft. 

Notice that in the formulation, which maximizes the time on the g-constraint , 
the optimization variable, T^, is also one of the independent variables. 

This is not uncommon, but is, nevertheless, confusing to new users, who are not 
familiar with optimization theory. 

As mentioned earlier, the complete SST0 entry problem is solved using the 
multiple run feature of POST. The first run is a one-dimensional iteration 
problem that determines the bank angle schedule required to achieve the proper 
pull-out conditions. The second run, initialized using the converged condi- 
tions from Run 1, simulates closed loop steering, first, along the heat rate 
boundary , and, second, along the acceleration boundary. When the time of de- 
parture is reached, the program uses the projected gradient algorithm to maxi- 
mize Tp, subject to the crossrange constraint. Notice that only two runs and 

one pass at the computer are required to efficiently solve this problem. This 
is because the external iteration (Run 1) and optimization (Run 2) occur at the 
"top" and "bottom" of the trajectory. The majority of the flight time is 
simulated only once using closed— loop steering to satisfy both environmental 
constraints . 

The data deck for this example is listed in Figure 32. 
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PSSEARCH 
SRCHM * 4, 

MAXITR = 3, 

NINOV * 1, 

INDVR * 5HDBANK, 

INDPH * 20., 

U * -71.1220821, 

PERT « .5, 

NDEPV * 1, 

DEPVR = 6HHEATRT, 

DEPPH * 40., 

DEPTL « .1, 

DEPVAL * 100., 
t 

PSGENDAT EVENTI1) * 1.0, 

7ITLE = 50HMAX HEAT RATE ^100. GAMMAI * -0.80, PERIGEE * 0 
TITLE C7) * 40HEAST MISSION . 

FESN * 100., 

NAXTIM * 4000. , 

ALTMIN * -110000., 

ALTMAX « 450000., 

NPC (II = 3, 

NPCC2I = 1, 

DT * 20., 

DTIMRtl > « 1., 

DTIMR (2 ) = 1., 

DTIMR (3 ) * 1., 

NPC 13 ) « 2, 

VELI = 25600., 

GAMMAI * -.8, 

AZVELI » 0., 

AZVELI * 90., 

NPC 1 4 J * 2, 

ALTITO = 400000., 

GCLAT >0., 

LONG * 0., 

J2 * 1.0827E-3, 

MU * 1 .4076469E+16, 

OMEGA = 7.292115E— 5 , 

PE * 20925722., 

RP * 20855560., 

NPC ( 5 J = 2, 0, 0, 2* 

LREF * 200. , 

RN *1 .0 , 

NPC <9 1 * 0, 

SREF * 12120., 

WGTSG * 523076., 

NPC (101*0, 

NPC (111 * 0, 

MONX(l) = 6HDYNP ,6HASMG ,6HHEATRT, 


Figure 32.- Input Data Deck for SSTO Sample Entry 
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MONY(l) * 6 HTIME , 6 HTIME , 6 HTIME * 

NPC (12) = Of 

IGUID * Of 0« 3* 

ALPHA = 30. f 

DALPHA * 30. » 

BNKANG * — 9C. , 

DBANK = -90. f 
NSPEC * If 

SPECI * 400000. t 

SPECK 6 ) * 100. * 

SPECI (7) * 310., 

NPC (15) = l.f 

PINC * 50., 

PRNT(91) » 6HTIMRF1 , 6HXNAX1 ,6HXMAX2 ,6HXHAX3 , 6 HTLPWT , 6 HUBAR , 
6 HT IHR F2 , 6 H YXMX 1 .6HYXMX2 ,6HYXHX3 ,6HTIHRF3, 6HSPECV1, 
6 HVELAD , 6 HH • 

* 

PSTBLMLT 

$ 

PSTAB TABLF * 6 HCDT , 1 , 6 HMACH , 9, 1,1,1, 


4.0, 

.430 , 

4.5, .420, 

5.0, 

.410, 

5.5, .400, 6.5* .390, 

8 . 0 , 

A 

.383, 

10 ., .375, 

30., 

.340, 

100., .340, 

» 

PSTAB TABLE = 

6 HCLT ,1 

96 HMACH 9 

8 , 1 , 1 , 1 , 

4.0, 

.622, 

^•7| • 6 00 9 

5.5, 

.585, 

6.5, .570, 7.5, .560, 

8.5, 

.550, 

10 • « »540 9 

30., 

.520, 



S 

PSTAB TABLF = 6HALPHAT, 1 ,6HMACH ,6, 1,1,1* 

0.,10.0, 1.0,10.0, 2.0,20.0, 4.0,20., 5.0,30.0, 100. ,30.0, 

ENDPHS *1, 

S 

PSGENDAT EVENT ( 1 ) * 10., 

CR1TR * 6HALT1T0, VALUE = 325000., 

NPC ( 1 ) * 0, 

TlHRF(l) * 0., 

TINRF (2 ) = 0., 

DT = 10., 

ENOPHS * 1* 

* 

PSGENDAT EVENT(l) * 20., 

CRITR * 6HTIMRF1 , 

VALUE * 0.0, 

T1MRF (1 ) * 0., 

MDL * 1 , 

DBANK * -71.1220821, 

DT * 2 . , 

ENDPHS * 1 • 

S 

PSGENDAT EVENT(l) =30., 

CRITR = 6HTIMRF1, 

VALUE = 10., 


Figure 32.- Continued 
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TIHRFCll * 0., 

IGUID * 0, 0, 0, 

HDL * 1 t 

ENDPHS * 1* 

* 

PAGENDAT EVENT * 35., 

CRITR = 6HVELAD * 

VALUE * -.1, 

ENDPHS * 1. 

A 

PAGENDAT EVENT = AO.. 1.* 

CRITR * 6HSPECV1. 

VALUE * 5.85. 

HDL * 1 , 

DT * 2.. 

TIHRF (1 ) * 0., 

IGUID ( 1 ) * 0. 

IGUID (2 ) * 1. 

IGUID (6 ) * 0, 

IGUIDC7) = 0, 

IGUIDC8) * A* 

KDG ( 3 I « *20. « 

KRG(3 ) = 300.. 

IDGFI3) = 1. 

DGF ( 3 ) = 6HHEATRT . 

TIHRF C3) * 0., 

A 

PATBLHLT 

$ 

PATAB TABLE = 6HGN0H3T . 1 .6HVELR . 7. 1.1. I. 

0., -70., 12000., -70,, 13500., -20., 

15500., -20., 18500., -A5., 2A000., -70., 30000., -70., 

A 

PATAB TABLE * 6HGNMX3T,0,0. , 

A 

PATAB TABLE = 6HGNMN3T, 0 ,-180. , 

A 

PATAB TABLE = 6HG0F3T ,0,100., 

ENDPHS =1, 

A 

PAGENDAT EVENT * 50., 0., 

CRITR c 6HASMG , 

VALUE = 1 • A , 

MDL s l , 

TIHRF (1 ) s o., 

TIHRF C3) * 0., 

KDG ( 3 ) s *200. t 

KRG ( 3 ) s 3000., 

IDGF(3) s 1, 

DGF (3 1 * 6HASHG , 

A 


Figure 32.- Continued 
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PATBLMLT 

F*TAB TABLE = 6HGN0M3T, 1 »6HVELR * 7,1,1, 1, 

0.» -70.* 12000.* -70., 13500., -20., 

15500. » -20. » 18500.* -45., 24000., -70..* 30000.* -70.* 

$ 

PJ TAB TABLE = 6HGNMX3T *0*0.* 

$ 

p$TAB TAeLE = 6HGNMN3T, 0 ,- 180 ., 

PSTAB TABLE = 6HG0F3T * 1 ,6HTIMRF 1*3* 1 * 1*1 * 

0., 1.4, 30., 1.4, 10.E10, 1.4, 

ENDPHS =1* 

$ 

PAGENOAT EVENT ( 1 1 * 60., 

CRITR * 6HTIMRF3, 

VALUE = 141., 

IGUID = 0,0,3, 

OALPHA =30., 

BNKANG = -70., 

OBANK = -26.7182816, 

T1MRF (3 1 = 0., 

MOL * 1 * 

OT =10., 

ENOPHS * 1, 

A 

PAGENOAT EVENT = 65., 1., 

VALUE = 90.0, 

CRITR = 6HAZVELR, 

CR1TR = 1HU, 

VALUE =0., 

IGUID (1 ) =0., 

1GUIC(2 ) = 1., 

IGUID(6 ) = 2., 

IGUI0(7» = 0., 

IGUID (8 ) = 1., 

BNKPC(l) = 0., 

ENOPHS =1, 

A 

PAGENOAT EVENT(l) =70., 

CRITR = 6HTIMRF3* 

VALUE = 20., 

MOL * 1 , 

TIMRF (3 ) = 0., 

OBANK * -23.6650460, 

ENOPHS = 1 , 

A 

PAGENOAT EVENT(l) * 75., 

CRITR = 6HTIMRF3, 

VALUE =20., 

MOL = 1, 


Figure 32.- Continued 
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TIMRF(3) 

* 0 • , 

DBANK » 

-11.1222061, 

ENDPHS 

* It 

% 


PSGENDAT 

EVENT * 76., 

CRITR 

= 6HTIHRF3 

VALUE 

= 20. , 

TIMRF (3 ) 

= 0 • , 

MDL 

* It 

DBANK * 

0.0, 

ENDPHS 

* It 

% 


PSGENDAT 

EVENT ( 1 ) = 

CRITR 

» 6HHACH 

VALUE 

= 5., 

MDL 

= It 

IGUID(l) 

=• 0. , 

IGUID(2 ) 

* 1 . t 

IGUIDC6I 

* 2. t 

IGUID(7 ) 

* 0., 

IGUID (8 ) 

* 1 . t 

BNKPC (1 ) 

* 0., 

DT * 5., 


S 


PiTBLMLT 



PiTAB TABLE * 

6HCDT 

,2 ,6HMACH , 

6HALPHA , 15 

t 14, 

1,1, 1,1 

,1,1. 

-90., 

0.00, .030, 

0* 20 ♦ 

♦ 0309 

0.60, 

.030, 0.70, 

.033, 

0.80, 

.039, 

0.90, .050, 

l.lOt 

.105 f 

1.20, 

.122, 1.30, 

.130, 

1.40, 

.125, 

1.60, .112, 

2*00 9 

.100 9 

5.00, 

.069, 10.0, 

.050, 

30.0, 

.026, 

o 

• 

o 

o 

0.00, .030, 

0*20 9 

.030, 

0.60, 

.030, 0.70, 

.033, 

0*80 9 

.039, 

0.90, .050, 

1.109 

.105, 

1.20, 

.122, 1.30, 

.130, 

1*40 9 

.125, 

1.60, .112, 

2.00 9 

.100, 

5.00, 

.069, 10.0, 

.050, 

30.0. 

.028, 

4*00, 

0.00, .042, 

0.20, 

.042, 

0.60, 

.042, 0.70, 

.045, 

0.80, 

.050, 

0.90, .067, 

1.10, 

.130, 

1.20, 

.142, 1.30, 

.145, 

1.40, 

.141, 

1.60, .126, 

2.00, 

.110, 

5.00, 

.070, 10.0, 

.050, 

30.0, 

.028, 

8.00, 

0.00, .055, 

0.20, 

.055, 

0.50, 
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Example 3. Orbital Maneuvers 


This example illustrates the application of POST to a finite-bum orbit 
transfer problem. The basic problem is to determine the optimal location, dur- 
ation, and attitude of two thrusting maneuvers that transfer an orbital vehicle 
(Transtage) from a near-Earth park-orbit to a geostationary orbit. Final weight 
is again used as the optimization variable. However, as demonstrated in Example 
1, this is the same as maximizing the deliverable payload. A simplified version 
of this problem could be easily set up using the impulsive option. The more 
detailed simulation is presented, however, to illustrate a number of specific 
flight phases that naturally occur in orbital analysis problems. These phases 
are: (1) non— Keplerian coast segments using the Krogh integrator, (2) attitude 

reorientation maneuvers that simulate the kinematics (not dynamics) of the RCS, 
(3) propellant settling via thrusting maneuvers, and (4) long-duration finite- 
thrust orbital maneuvers. 

The basic transfer geometry is illustrated in Figure 33. As indicated, the 
problem starts at booster burnout, which occurs 483.893 second into the flight. 
The first phase. Event 1.0 to Event 125.0, is a non— Keplerian coast to the 
first equatorial crossing. The occurrence of the equatorial crossing is de- 
mined by specifying "latitude equals zero" as an event criteria. This is accom- 
plished via the $GENDAT input shown below. 

P$ GEN DAT 
EVENT = 125. , 

CRITR = 5HGCLAT , 

VALUE =0.0, 

ENDPHS = 1, 

$ 

Event 125.0 is used to "trigger" a sequence of (1) reorientation, (2) pro- 
pellant settling, and (3) main engine start-up transient phases. The primary 
main engine bum starts at Event 148.0 and ends some 325 seconds later at Event 
150.0. Event 150.0 also initiates the coast -to-apogee defined by Event 260.0. 

At apogee the same basic sequence of maneuvers are used again to simulate the 
final circularization maneuver, which ends at Event 280.0. A long coast tra- 
jectory is then propagated until the vehicle reaches its fourth equatorial 
crossing at the final event. Event 300.0. 

The POST input deck for this problem is shown in Figure 34. Comments are 
placed on each important data card to assist the reader in understanding the 
input set up. The ability to place comments on actual data cards is a feature 
unique to the POST NAMELIST processor, and is not available on most standard 
system NAMELIST routines. 

The targeting and optimization formulation and input set up used in this 
case deserves some explanation. First, the optimization variable is WEIGHT at 
Event 300.0. This is also the WEIGHT at Event 280.0 because the variable WEIGHT 
does not change during the coast phase between Events 280.0 and 300.0. New 
users often think this practice to be inefficient. However, in this case it is 
not because the orbital conditions (constraints) are also defined at the fourth 
equatorial crossing, and so the trajectory must be propagated to Event 300.0 to 
calculate the dependent variables. 
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| represents the crossing condition 



Event 300.0 - End of Simulation at 4th Equatorial Crossing 
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Figure 33.- Orbital Transfer Geometry 
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Second, there are five orbital "target" conditions defined that describe 
the desired final orbit. The first four constraints are the standard Keplerian 
orbital elements specifying semimajor axis, eccentricity, inclination, and 
argument of perigee. The last constraint is longitude at Event 300.0. Now 
Event 300.0 defines the fourth equatorial crossing, which in this case is an 
ascending crossing. As a result, the final constraint is mathematically equal 
to the longitude of the ascending node. At first this seems a rather obtruse 
procedure for computing the longitude of node. However, this special technique 
is generally required to eliminate the problems encountered in targeting to 
very low inclination orbits, where as a result of small changes in the maneuver, 
the ascending node can discontinuously change by 180 deg. This is because 
during a long thrust maneuver a "seemingly" ascending node can be changed to a 
descending node with only minor changes in the yaw-turn steering parameters. 

Third, there is a total of nine independent variables defined in this sample 
case. The first five of these variables are, because of their event numbers, 
associated with the first main engine apogee-raising bum. The first control 
variable in the event criteria at Event 140. Review of the data deck reveals 
that CRITR at Event 140.0 is latitude. Thus, this first control variable is the 
latitude at which the apogee raising maneuver is started. The second and third 
control variables are the initial relative pitch and yaw attitude angles, re- 
spectively. The fourth variable is the pitch rate during the burn, and the 
fifth variable can be shown to be the bum time of the main engine. The final 
four control variables are similarly defined, and apply to the final circular- 
ization maneuver. Notice in the final maneuver a constant relative attitude is 
maintained because the pitch rate coefficient PITPC2 has been omitted from the 
final set of control parameters. 

The mathematical formulation of this orbital transfer problem may be sum- 
marized as follows: 

Determine the values of the nine control parameters 

u = 6 l» ♦!» 6l» t b1 » <f>2> e 2» ^2» T b2 ) (41) 


that maximize: 

subject to the final orbital conditions (constraints); 


M*f) 

- 13 893 648 ft 

e (* f ) 

- 0.030 = 0, 

H' 

rr 

■“h 

- 5.0 deg = 0, 

u 

3 

- 180 deg = 0, 

/t , i 
R \ f J 

| - 58.96 deg = 0 
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where $ is geocentric latitude, 0 is relative Euler pitch angle, <l> is relative 
Euler yaw angle measure from the azimuth of the inertial velocity vector, 9 
is the time derivative of 0, T„ is the bum time of the main engine, a is semi- 

major axis, e is eccentricity, i is inclination, u> is argument of perigee, 0 R 

is the longitude of ascending node, and W is total spacecraft weight. The 
subscripts 1, 2, and f denote the first and second maneuvers and the final con- 
ditions, respectively. The relationship between this mathematical formulation 
and the required $ SEARCH inputs can be distilled by studing Figure 34. There 
are a number of additional inputs in $ SEARCH that are optimization algorithm 
control parameters, and are not related to the problem formulation. Brief de- 
scriptions of these parameters are given in the comments on the individual data 
cards. More complete definitions are contained in reference 1. 


Example 4. Hypersonic Aircraft Point Performance 

This example is presented to illustrate the application of the POST pro- 
gram to aircraft point performance problems. This is a unique application of 
POST in that the trajectory time history is not generated. Instead, the com- 
plete optimization takes place at a single fixed point in time— hence, the 
term point performance. See Figure 35. 

The basic problem is to determine the maximum cruise velocity of a hyper- 
sonic aircraft such as the X-24C. The maximum velocity, of course, depends on 
the configuration and the particular propulsion system used. For a particular 
configuration and engine, the problem may be stated more precisely as follows: 

Determine the cruise altitude, velocity, and angle of attack h, V R , a 
to 

maximize: 

subject to: 


The constraints on the total derivative of relative velocity and flight path 
angle are used instead of the standard cruise conditions of "thrust equals 
drag' 1 and "lift equals weight." This is due to two important accelerations 
that cannot be ignored in the hypersonic flight regime: (1) the thrust com- 

ponent in the lift direction, and (2) the centripetal acceleration. 

The complete input deck for this example is presented in Figure 36. There 
are only two events in this case: ' EVENT = 1.0, the first event, and EVENT = 
100.0, the final event. The final event criteria, CRITR = 5HTDURP , is the 
elapsed time between EVENT - 1.0 and EVENT = 100.0. Notice that the value of 
the final event criteria is zero, i.e., VALUE = 0.0. This means that zero time 
evolves between events 1.0 and 100.0. This type of event is referred to as a 
"zero-length event" in POST terminology. 
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Figure 35.- Mission Profile and Vehicle Sketch for Hypersonic Point Performance 
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Figure 36.- Input Data Deck for Hypersonic Aircraft Sample Case 
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The iterative optimization begins by the user estimating an initial "guess" 
for the independent variables. 

u = (h, V R , a) = (90 000., 4500., 10.). (42) 

This initial guess, input in $SEARCH (see Figure 36) generally does not satisfy 
the constraints (cruise conditions). The projected gradient algorithm then 
determines the proper correction to this initial guess to satisfy the constraints. 
This correction is based upon numerical approximations to the partial derivatives 
of the constraints. This set of partials, called the sensitivity (or more pre- 
cisely the Jacobian) matrix is given as 


fsi Hv ir! 

l J2x3 9(h, V R , a) 


After a feasible solution is attained, that is, one that satisfies the constraints, 
the algorithm is designed to "move along" the constraint manifold and improve the 
cruise speed at each iteration. Conceptually, this is accomplished by projecting 
the gradient of cruise speed to the plane defined by the normals to the con- 
straints. This mathematical operation determines the direction of "steepest 
ascent (descent)" in a plane that approximates the constraints. However, be- 
cause of nonlinearities in the constraints, as the control parameters are varied 
along this direction of search the cruise condition becomes increasingly violated. 
The extent to which the control parameters are incremented is then determined by 
maximizing a composite function (called estimated net performance), which is the 
sum of the cruise speed and an estimate of the loss in cruise speed associated 
with this violation of the cruise conditions. After each such optimization step 
a constraint restoration step is taken to remove this induced constraint error. 

As a result, once a feasible solution is found, the optimization proceeds in a 
sequence of two steps: first, optimization; and second, constraint restoration. 

Significant, however, is the fact that in most normal restoration steps, the 
sensitivity matrix is not recomputed. This procedure substantially reduces the 
overall computer solution time required to optimize the problem. A brief sum- 
mary of this iteration process is given in Table 14. Notice that after a feas- 
ible solution is achieved in iteration 1, all subsequent iterations satisfy the 
cruise equations. 

This simple application of POST solved this "point" performance problem 
in a matter of only a few seconds on the CDC 6500. As a result, it eliminates 
the requirement of plotting families of thrust required and thrust available 
curves to determine graphically the cruise performance of a high performance 
aircraft configuration. 
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Example 5. Trajectory Decomposition 


The optimization of an eastern hemisphere geostationary mission using the 
decomposition option is presented as Example 5 • This mission, illustrated in 
Figure 37, has two basic flight regimes and is typical of the problems that can 
be efficiently optimized with the decomposition approach. It is important to 
note that the decomposition option is contained in a special version of 3D POST, 
However, the utilization procedures are available in reference 1, A general 
problem statement and the two formulations are given in Table 15. The standard 
nondecomposition formulation is a constrained discrete parameter optimization 
problem containing eight constraints and ten independent variables. Let this 
problem be denoted as PI (8x10), The data deck is presented in Figure 38. The 
Inputs required by the decomposition appear first in the deck. The remainder 
of the data deck is the same as the standard nondecomposition input. As indi- 
cated, the complete mission was decomposed to three trajectory segments: 

Ascent; 

Apogee Raising Bum; 

Final Circularization and Plane Change. 

This problem could be decomposed differently; for example, the two orbital seg- 
ments could easily be combined in a single subproblem. The master problem, as 
formulated, has only one constraint and three control variables; let the prob- 
lem be denoted as MPl(lx3). It is important to note that in MPl(lx3) two of 
the three control variables are constraint values in the subproblems. This en- 
ables the master algorithm to optimize the burnout velocity of the ascent leg 
and the inclination of the transfer orbit. These are two important geosta- 
tionary mission design variables. The constraint that zero propellant remains 
at final burnout is really the active inequality constraint W pr >_0. This is 

to ensure that the program will not simulate consumption of more propellant 
than can be loaded in the tanks. This constraint is defined at the master level 
because it couples the subproblems due to the fact that Stage III is used in 
all three of the trajectory segments. The subproblems, represented as SPl(3x3) , 
SP2(2x2) , and SP3(4x4), are all relatively small targeting problems. Clearly, 
this sequence of subproblems is much easier to solve than PI (8x10) • However, 
a computational tradeoff is involved because the sequence of subproblem 

/SP : i * 1. 2, 3\ must be solved on every trial step in univariate searches 
\ i / 

used to solve MPl(lx3). This fact emphasizes the importance of rapid subprob- 
lem solution. 

Computational results for each formulation are presented in Table 16 and 
17. The accelerated gradient projection option was used to solve the standard 
formulation. On iterations 3 and 7, it had to be restarted and did not con- 
verge in ten iterations. In comparison, the decomposition algorithm did not 
require manual restarting and converged in only six iterations. The increased 
reliability is probably more significant than the reduction in the number of 
iterations. This is because each iteration at the master level represents con- 
siderably more computational effort than does each iteration in the standard 
formulation. 
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(l) Stage II Shutdown 


( 5 ) Stage III First Bum Start 

(D Stage III First Burn End 



Typical Eastern Hemisphere Geostationary Mission in Earth Centered 
Relative Coordinates 


Figure 37.- Mission Profile for Decomposition Sample Case 


Example 6. 6D0F Space Shuttle Entry Simulation 

A typical data deck for a 6D0F closed-loop simulation of the Space Shuttle 
entry trajectory is presented. The changes, additions and modifications re- 
quired to convert a "typical" 3D POST Space Shuttle entry data deck to this 6D 
POST data deck are highlighted in Figure 39. As indicated, the primary changes 
are related to the additional data cards required for (1) mass properties, (2) 
aerodynamics, and (3) the flight control system (FCS) . More significant than 
these minor data additions are the vehicle-dependent subroutines that must be 
supplied by the user to model the FCS. Development, coding, and checkout of 
these subroutines represent the vast majority of the engineering effort in using 
6D POST. However, this requirement is certainly not unique to 6D POST, but 
rather is required by all 6D0F simulation programs. The key issue then is 
really related to the ease in which these routines can be implemented into the 
program. These procedures, described in reference 5, have been made as simple 
as possible to reduce the "start-up" costs of using 6D POST. This simplistic 
simulation of Space Shuttle entry requires about 100 000 octal cells of com- 
puter memory, and executes in about three-times realtime on the CDC 6500. 
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rmine the attitude steering history and burn tiroes that maximize the payload a Titan IIIC class vehicle can deliver 
n eastern hemisphere geostationary orbit. Assume that the elements for the optimum park orbit and transfer orbit 
unknown and are to be determined. 
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TABLE 16.- NONDECOMPOSITION ITERATION SUMMARY 
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COMPUTER REQUIREMENTS 


Machine Configurations 

POST was originally developed for the CDC 6000 series computers. 3D POST 
has since been converted to the IBM 370/165/OS and the UNIVAC 1108/EXEC 8 com- 
puters. Both 3D POST and 6D POST are coded exclusively in FORTRAN IV and are 
compatible with the FORTRAN IV EXTENDED compiler. The 6D POST is not available 
on any computer /system other than the CDC 6600/NOS. The minimal computer re- 
quirements for 3D and 6D POST are given in Table V-l. 


TABLE 18.- MINIMAL POST COMPUTER REQUIREMENTS 


Computer 

Operating 
1 System 

Precision 

3D Storage, 
Words 

6D Storage, 
Words 

CDC 

6400 

6500- 

6600 

NOS, SCOPE 
3.4 

Single 

140 000 * 

O 

140 000 o 

o 

IBM 

370/165 

OS 

Double 

50 000 

70 000 

UNIVAC 

EXEC 8 

Double 

55 000 

65 000 

*Nonoverlayed 

version. 





As with any larger computer program, special operating system-dependent 
techniques can be used to reduce core requirements or increase execution speed. 
For example, the overlay structure can be tailored to any given computer/system 
to reduce core size and execution time. The specific compiler used can also 
significantly impact execution speed. For example, the CDC 6600 POST executes 
about two times faster using the FTN compiler than it does using the RUN com- 
piler. These system-dependent program tradeoffs and modifications can easily 
be made once the standard version is operating correctly on any particular com- 
puter system. 


Computer Precision 

There are two important numerical techniques used in POST that require con- 
siderable computational precision to function properly. These techniques are 
(1) numerical differentiation, and (2) numerical integration. The most diffi- 
cult computational problem encountered in POST combines both numerical tech- 
niques, that is, the approximation of trajectory partial derivatives by numeri- 
cal differentiation of numerically integrated trajectories. This numerical 
process dictates the computer precision requirements for POST. Computational 
experience indicates that a 64-bit word is usually adequate for the numerical 
differentiation of various trajectory variables. This means that double 
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precision is required to run POST on any computer that has a 32-bit single pre- 
cision word. As a result, the program must be executed in double precision on 
IBM and UNIVAC computers. Single precision is adequate on all CDC 6000 series 
computers because CDC has a 60-bit single precision word. 

To numerically differentiate the dependent variables (target conditions, 
inequality constraints, and the optimization variable) with respect to the in- 
dependent variables, the user must input the perturbation to be used for each 
independent variable. The user must also select the integration method and the 
initial integration step size. Generally, the numerical truncation error asso- 
ciated with the numerical integration process tends to cancel itself in the dif- 
ference quotient; however, propagation of local round-off error does not cancel. 
Fortunately in most trajectory problems, the number of integration steps is 
sufficiently small so round-off error is not a major problem. As a result, the 
selection of the integration method and step size, although important in terms 
of run time and accuracy, is not critical to the numerical differencing tech- 
niques employed. Critical, however, is the proper selection of the control vari- 
able perturbations that are used in the divided difference formula. If the 
perturbations are excessively large, accuracy is lost because of truncation 
error in the first and second difference formulas. If the perturbations are too 
small, accuracy is lost due to subtraction of the nominal value from the per- 
turbed value on a finite word length computer. Problems associated with the 
word length can, and have for the most part, been eliminated by using either CDC 
computers or by using double precision arithmetic on computers with smaller 
single precision word length. However, control of the truncation error in the 
difference formula is difficult and numerical experience on a given problem/ 
computer/operating system is sometimes required. To alleviate this problem an 
automatic perturbation step size controller is incorporated in the latest ver- 
sions of the programs. This routine requires an initial guess for the perturba- 
tions in each independent variable. The initial guess is used unless the re- 
sulting divided difference quotient is judged to be sufficiently inaccurate as 
to require modification. Modification is accomplished by first rerunning the 
perturbation with its sign changed to obtain a symmetric difference for those 
variables with incorrect perturbations. The symmetric differences are then 
used on the present iteration. On the next iteration, the perturbations are 
adjusted to ensure the proper amount of variation in the dependent variables 
to secure adequate accuracy in the sensitivities. 


Runtime 

The runtime, as measured by the Central Processor Unit (CPU) clock, is a 
key input to all computer charging algorithms. As a result, it is important to 
understand what factors contribute to the CPU time requirements when running 
POST. Understanding will enable users to make reasonable estimates of computer 
budgets associated with using POST programs. 

The CPU time for a POST run depends on numerous factors. The most impor- 
tant are as follows: 
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Computer hardware and system software characteristics; 

Type of trajectory and vehicle being simulated; 

Accuracy requirements; 

Number of degrees of freedom in the targeting/optimization formulation. 

The computer and the operating system used clearly play a primary role in 
determining the runtime of any program. In addition, the CPU time is not the 
only variable used in determining the cost-effectiveness of a given program/ 
computer system interface. For example, the CPU time required to run a POST 
input deck on the CDC 6600 /NOS /FTN may be only a fraction of that required to 

run the same deck on the CDC 6400/SCOPE 3.4/RUN; however, it may be less ex- 

pensive to use the older and slower computer because of the charging algorithms 
involved and the demand on each system. The point being made is that the only 
practical way to estimate computer costs is through computational experimenta- 
tion with POST on your computer system/problem. Generally, minor changes can 
be made to the program executive structure (overlays) to eliminate any severe 
incompatibilities between POST and the host computer system. 

A second key factor impacting runtime is the type of trajectory and ve- 
hicle being simulated in any given problem. For example, exoatmospher ic tra- 

jectories can be integrated more rapidly than atmospheric ascent of entry 
trajectories. Also problems that require large amounts of tabular data to de- 
scribe the mass properties, the aerodynamics, and the propulsion system will 
require longer corresponding runtimes. The total flight time of the trajectory 
being simulated is a prime driver in the computational costs. Experience indi- 
cates that long duration flights generally require more computer time than 
shorter flights because of the direct increase in the number of required inte- 
gration steps. The amount and frequency of printout data requested are also an 
important consideration that can always be controlled by the user. In extreme 
cases, runtime can double or even triple when the full 198 variable printblock 
is requested at every integration step. Thus, to reduce runtime, only the mini- 
mal amount of data required to accurately model the vehicle should be input and 
only the minimal amount of data required to interpret the results should be 
output . 


Accuracy, which is directly related to the selection of the integration 
method and the integration stepsize, is an important consideration in every 
problem. Generally speaking, the standard fourth order Runge-Kutta algorithm 
with the appropriate stepsize is the most cost effective numerical integration 
method; however, under special circumstances higher order methods may be faster 
For example, Keplerian orbits can be integrated ten times faster with the Krogh 
method than with Runge-Kutta. The best approach to determining the impact 
stepsize on accuracy and CPU time is to run several single pass trajectories. 

A compromise stepsize can then be determined from plots of integration error 
and CPU time versus stepsize. 
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The final factor influencing runtime is the number of degrees of freedom 
difference^" 8 and optimization formulation. The number is computed as the 
n.™w f b ® tW6en the number of ^dependent variables, m, and the average 
number of active constraints, n^ The number of degrees of freedom can be use, 

to establish an estimate of the number of iterations, N, required by the accc) 
erated projected gradient algorithm to achieve the optimum trajectory. This 
estimate is given by the formula J y inis 


N = m - n + K 
a 


( 1 / 

where K is the number of iterations required to achieve a feasible, targeted 

condi t-i ° ry * a* eXa " iple ’ if there a re ten independent variables, three target 
conditions, and two inequalities, one of which is active, then 8 


N = (10 - (3 + 1)) + k = 


6 + K 


( 2 ) 


to ?r rally 3 reas ° nable estimate of the total number of iterations required 
to achieve an optimum trajectory. The number of iterations required to achiev 

J, trajector y’ K ’ var ies ** a function of the initial guess. Jo! T 
_easonably accurate initial guess, K can be estimated as the integer part of 
n a /2. In the case where there are no inequality constraints, then K is equal 

to the number of target conditions defined by the user. 3 

The total CPU time can then be approximated as 

CPU = N(T/I) 

where T/l is the average time required to make a single iteration 
estimated in terms of the single pass CPU time, x, as 

T/I = (1 + m + 6) x 


( 3 ; 

T/I can be 


where 1 + m trajectories are required to obtain the sensitivity matrix, and si 
rajectories are required (on the average) to perform the univariant searches 
(minimum is two and the maximum is 10) . This estimate for CPU per iteration i- 
conservative due to the fact that POST does not integrate the complete trajec- 

cnll ,° m perturbed runs. As a result, it is generally more ac- 

curate to determine T/I directly from actual computations. For this reason, tl 
CPU time used per iteration is computed and included in the iteration summary 

aS r f r /ITR - calculation of CP/ITE retires a machina-dapondent 

subroutine that usually must be modified for each computer/system. 
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